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1. Background 

The Mars Global Reference Atmospheric Model (Mars-GRAM) was developed 
(Johnson et al., 1989; Justus 1990, 1991) as an engineering-oriented, empirical model of the 
Mars atmosphere. A complete history and description of the model, through version 3.34, was 
recently given in the Mars-GRAM programmer's guide (Justus et a]., 1996). Mars-GRAM is 
based on surface and atmospheric temperature data observed during the Mariner (orbiter) and 
Viking (orbiter and lander) missions and on surface pressure data observed by the Viking 
landers. At the higher altitudes (above about 120 km), Mars-GRAM (through version 3.34) 
was based on the Stewart (1987) model for the global-mean thermosphere. 

Mars-GRAM provides both mean and mountain-wave perturbed atmospheric density 
for any location (height, latitude, and longitude) and time (seasonal and diurnal). Other 
atmospheric variables provided include atmospheric temperature, pressure, and wind 
components. Dust storm effects are included for all atmospheric parameters controlled by 
user-selected options. The model also includes the option to simulate either local-scale or 
global-scale dust storms and density perturbations from tidal waves, computed by the Zurek 
wave model, Pitts et al. (1990?). Other recent features include (1) a limitation, based on 
atmospheric stability considerations, for the magnitude of the mountain-wave density 
perturbations, (2) comparisons of density, temperature, and pressure with the COSPAR 
reference atmosphere (Pitts et al., 1990?), and (3) a new method for estimating the diurnal 
range of the surface temperature based on the diurnal variability of the surface-absorbed solar 
energy. Technical descriptions of all these features are given in the programmer's guide 
(Justus et al., 1996). 

The original Stewart (1987) thermospheric model, designed to estimate global mean 
conditions only, has two major shortcomings for use as the thermospheric portion of Mars- 
GRAM: (1) does not produce realistic variations of temperature and density with latitude and 
time of day (or longitude), (2) is not capable of providing realistic horizontal gradients of 
pressure, from which thermospheric winds are estimated. The additions to Mars-GRAM 
discussed here are designed to overcome these deficiencies. 

Mars-GRAM was developed from parameterizations to atmospheric data observed by 
the Mariner (orbiter) and Viking (orbiter and lander) missions and thus, is representative of 
the Mars atmosphere during the 1970's. Recent observations by Clancy et al. (1990) indicate 
that, in response to atmospheric cooling as the dusty atmosphere has cleared in the last two 
decades, current Mars temperature profiles are distinctly cooler than those observed in the 
Viking era. Comparisons of Mars-GRAM 3.34 mid-latitude average temperature profiles with 
recent data from Clancy (provided by Rich Zurek, private communication), indicate about 20 
K cooling in the 40-50 km height range, about 15 K at 20-30 km, but little change in the 5-10 
km region. This temperature change could have significant effects on atmospheric density at 
high altitudes and, as pointed out by Clancy et al. (1990), is quite important in planning for 
Mars missions that involve aerobraking. 

To accommodate these possible effects of climate shifts, Mars-GRAM version 3.4 
allows the user to select climate modification factors for the temperature profile between the 
surface and 75 km. Additionally, the user may adjust at run time relevant values of the global 



mean thermosphere (exospheric temperature, and temperature and height of the base of the 
thermosphere). 

Proper computations of variation with time of day depend on accurate estimates of 
longitude of the sub-solar point on the surface of Mars. In the original (1989) development of 
Mars-GRAM an empirical fit for sub-solar longitude versus time was derived from Mars 
ephemeris data for 1984 through 1988. This empirical fit was incorporated into the "ORBIT" 
subroutine. The original ORBIT subroutine used a time variation for sub-solar longitude 
based on a Fourier harmonic series with a basic period of 696 days. For Mars-GRAM version 
3.4, we obtained (Myles Standish, Jet Propulsion Laboratory, private communication) a set of 
Mars ephemeris data from 1984 through 2003. A Fourier fit with this longer time series 
revealed that the sub-solar longitude was reproduced significantly better with two harmonic 
series, one with the basic period of the Mars orbit (687 days) and another with 777 days. This 
new, dual-period parameterization was incorporated into the ORBIT subroutine of Mars- 
GRAM 3.4. 
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2. Methods for Thermospheric Calculations 


2.1 Overview of Thermospheric Model Methods 


Five basic parameters are used to prescribe the height variation of pressure, density, 
and temperature in the Stewart (1987) thermospheric model: (1) pressure (PRESSF) at base 
of the thermosphere, assumed to be the fixed value 1.26 nanobars (1.26 x 10' 4 N/m 2 ), (2) 
height (ZF) of base of the thermosphere, (3) temperature (TF) of base of the thermosphere, (4) 
exospheric temperature (TINF), and (5) a parameter (SCALE) that describes the height 
variation of temperature between base of the thermosphere and the exosphere. In the original 
Stewart thermospheric model, ZF, TF, TINF and SCALE are considered global average values 
dependent on either the 10.7-cm solar flux (F10.7) or the heliocentric distance of Mars from 
the Sun, in Astronomical Units (RAU). In the new thermospheric model, these parameters are 
also functions of latitude and time of day (longitude). Details of computation of the height 
variation of temperature, pressure, and density are identical to the original Stewart 
thermospheric model. 

For an altitude (Z) above the base of the thermosphere of ZZF = Z - ZF, located where 
the radius of the Mars reference ellipsoid is RO (and RF = RO + ZF), the vertical temperature 
variation is computed by 


YSC = ZZF * RF / (RF + ZZF) 

TZ = TINF - (TINF - TF) * EXP(-YSC / SCALE) 


THRM108 

THRM109 


(Letters and number at the right are line numbers of the FORTRAN source code of the 
Mars-GRAM program.) Scale height (HH) and partial pressure (PRZ) for each 
thermospheric constituent are computed by 

HH(I) = BK * TINF / (GF * DM(I)) / 1 . 0E5 THRM111 

PRZ(I) = PRESSF*FF ( I ) *EXP (-YSC/HH ( I ) - (SCALE/HH(I) ) *ALOG (TZ/TF) ) THRM113 

where FF is a fractional composition for each constituent. This fairly simple analytical 
expression for pressure variation can be used because the particular functional form 
assumed for temperature variation (code line THRM109) yields a closed-form solution for 
the vertical integral of the hydrostatic equation. Note that the decrease of gravity with 
height (inversely proportional to the square of radial distance) is explicitly accounted for in 
this hydrostatic integral solution (code line THRM1 13). 

Density is computed from the total pressure (sum of the partial pressures) and 
temperature, using the perfect gas law relation. The problem of specifying temperature, 
pressure, and density at any height, latitude, time of day, time of year, and solar activity 
condition, becomes specifying dependence of ZF, TF, TINF, and SCALE on these factors. 
2.2 Dependence on Heliocentric Distance and Solar Activity 
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A collection of thermospheric parameter estimates, observations and model results, is 
assembled in table 2.1. The data were collected and given in model results or measurement- 
model comparisons by Bougher et al. (1988, 1990, 1993) and Bougher (1995). All parameter 
values in table 2.1 are global mean conditions. 


Table 2.1. Global mean values of thermospheric parameters from various measurements or 
model results [Bougher et al. (1988, 1990, 1993) and Bougher (1995)]. R is the heliocentric 
distance of Mars, TINF is exospheric temperature, TF and ZF are (respectively) the 
temperature and height of the base of the thermosphere, and SCALE is the height variation 
parameter for thermospheric temperature. 


Data Source 

Mariner 4 
Mariner 6-7 
Mariner 9 (primary) 
Mariner 9 (extended) 
Viking Lander 1 
Viking Lander 2 
Phobos (DMED) 

Phobos (DMAX) 

Phobos (DMIN) 

MTGCM, MGS97L 
MTGCM, MGS98L 
MTGCM, MANCOO 
MTGCM, MGS97T 


10.7-cm Solar Flux 


R (AU) 

1 AU 

Mars 

1.55 

77 

32 

1.45 

206 

98 

1.48 

110 

50 

1.62 

115 

44 

1.62 

68 

26 

1.63 

74 

28 

1.53 

150 

64 

1.67 

150 

54 

1.38 

150 

79 

1.38 

130 

68 

1.66 

130 

47 

1.53 

150 

64 

1.40 

80 

41 


TINF 

TF 

ZF 

SCALE 

(K) 

(K) 

(km) 

(km) 



117 


240 

170 

129 


220 

160 

123 




115 


180 

125 

125 


180 

125 

121 


215 


125 


200 


115 


230 


135 


228 

147 

128 

15.0 

196 

141 

115 

12.9 

229 

147 

122 

14.6 

185 

139 

127 

12.4 


Regression relations versus both heliocentric distance (R) and 10.7 cm solar flux at 
Mars (FI 0.7) were tested. The resulting best fit relationships to describe the variation of the 
global mean thermospheric parameters are the following: 

<TINF> = 156.3 + 0.9427 F10.7 

(2.1) 

<ZF> = 197.94 - 49.058 R 

(2.2) 

<TF> = 113.7 + 0.5791 F10.7 

(2.3) 

<SCALE> = 8.38 + 0.09725 F10.7 

(2.4) 


where the angle brackets denote global average values. 
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Figures 2.1 through 2.4 illustrate the regression relations ( 2 . 1 ) through (2.4) versus 
data from table 2. 1 . The example point (labeled VL-1) in each figure is the data point for 
the Viking 1 Lander (VL-1) case. The figures show the VL-1 point is reproduced fairly 
well for <TINF> and <TF> regressions, but the VL-1 point is a significant "outlier" in the 
<ZF> regression. Thus, in order for the model to reproduce specific cases, it may be 
necessary for the user to adjust the regression values by appropriate amounts (e.g., an 
increase of about 8 km in <ZF> to replicate the VL-1 case). Note that the VL -1 case was 
also a significant outlier in the regression plot Stewart (1987) derived for ZF versus 
heliocentric distance (from which he derived a 1/R regression relation for ZF). The user- 
input adjustment parameters (deltaTEX for exospheric temperature, and deltaTF and 
deltaZF for temperature and height of the base of the thermosphere) are discussed more 
fully below and in Appendix A. 

The regression relation for global mean exospheric temperature [equation (2.1) and 
figure 2.1] gives exospheric temperatures significantly lower than the original Stewart 
regression. However, the new exospheric temperature regression is quite similar to the 
curve described by the lower boundary of the shaded portion of figure 9 of Bougher et al. 
(1990), which also applies to global mean exospheric temperature. [The middle and upper 
curves within the shaded area of Bougher’s figure 9 are for mid-aftemoon and dayside- 
mean exospheric temperatures, not global mean values.] Note that equation (2.3) implies 
a dependence of T F on solar activity (F10.7), rather than an orbital radius (R), as assumed by 
Stewart. Residual error from regression via equation (2.3) is 8.1 K (essentially the same as 
residual error from multiple regression against both F 10.7 and R). Residual error from 
regression against R alone (yielding <T F > = 267.13 - 80.905 R) is 13.9 K. The physical 
reality of the dependence of T F on F10.7, implied by equation (2.3), is still being evaluated. 

2.3 Dependence on Latitude, Sub-Solar Latitude, and Time of Day 

Four sets of output from the Mars Thermospheric Global Circulation Model 
(MTGCM, Bouger et al., 1990) were used to derive the dependence of the thermospheric 
parameters on latitude, sub-solar latitude, and time of day. Parameterization relations for 
the various thermospheric parameters were derived and incorporated into a new subroutine 
(Thermpar, line numbers denoted TPAR, see Appendix A). The four MTGCM cases used 
were MGS97L , MGS98L", "MANC00" and "MGS97T"; characteristics are given in 
table 2 . 1 . 

In the subroutine Thermpar, the local (latitude and time of day dependent) value of 
exospheric temperature (TINFO, K) is computed from the global average value (<TINF> = 
Tbar) by the program steps: 


C... Zonal average exospheric temperature (K) versus latitude 
Tavg = TbarMl. + 1 . 369E-4*sunlat*lat) 

C... Phase angles (hours) for local solar time variation 
tl = 13.2 - 0 . 00119*sunlat* lat 
t2 = 9.4 - 0 . 0023 l*sunlat* lat 
C... Amplitude factor for local solar time variation 

cphi = Cos (pil80* (lat + sunlat)/(l. + LATMAX/90.)) 

C... Exospheric temperature (K) versus local solar time 

TINFO = Tavg* ( 1 . + 0 . 22*cphi*Cos <pil80*15 . 0* (LST-tl) ) + 

Sc 0 . 04*cphi*Cos (pil80*30 . 0* (LST-t2 ) ) ) 


TPAR 

47 

TPAR 

48 

TPAR 

49 

TPAR 

50 

TPAR 

51 

TPAR 

52 

TPAR 

53 

TPAR 

54 

TPAR 

55 

TPAR 

56 
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where "lat" is the local latitude (in degrees), "sunlat" is the latitude (in degrees) of the sub- 
solar point on the Mars surface, and "LST" is the local Mars solar time (in Mars hours = 
l/24th sol). 

Local height of the base of the thermosphere (ZFO in km) is computed from the 
global average value (<ZF> = Zbar) by program steps: 


c. . . 

Latitude variation factor 

TPAR 

59 


factlat = ( sunlat /LATMAX) *{lat/77.5)**3 

TPAR 

60 

C. . . 

Zonal average base height (km) versus latitude 

TPAR 

61 


Zavg = Zbar + 4.3*factlat 

TPAR 

62 

c. . . 

Amplitudes for local solar time variation 

TPAR 

63 


Al = 1.5 - Cos (pil80*4 . 0*lat) 

TPAR 

64 


A2 = 2 . 3* (Cos (pil80* (lat + 0 . 5*sunlat) ) ) **3 

TPAR 

65 

c. . . 

Phase angles (hours) for local solar time variation 

TPAR 

66 


tl = 16.2 - (sunlat/ LATMAX) *Atan (pil80*10 . 0*lat) 

TPAR 

67 


t2 = 11.5 

TPAR 

68 

c. . . 

Base height of thermosphere (km) versus local solar time 

TPAR 

69 


ZFO = Zavg + Al*Cos (pil80*15 . 0* (LST-tl) ) + 

TPAR 

70 


& A2*Cos (pi!80*30 . 0* (LST-t2) ) 

TPAR 

71 


Local temperature (TFO, K) at the base height of the thermosphere is computed 
from its global average value (<TF> = Tbar), by subroutine steps: 


C... Zonal average temperature at thermosphere base (K) vs. latitude TPAR 74 


Tavg = TbarMl. + 0 . 186*factlat) TPAR 75 

C. . . Amplitudes for local solar time variation TPAR 76 

Al = 0.06 - 0 . 05*Cos {pil80*4 . 0*lat) TPAR 77 

A2 = 0 . 1* (Cos (pil80* (lat + 0 . 5*sunlat) ) ) **3 TPAR 78 

C... Phase angles (hours) for local solar time variation TPAR 79 

tl = 17.5 - 2.5* ( sunlat /LATMAX) *Atan(pil80*10 . 0*lat) TPAR 80 

t2 = 10.0 + 2.0* (lat/ 77 . 5) **2 TPAR 81 

C... Thermosphere base temperature (K) versus local solar time TPAR 82 

TFO = Tavg* (1.0 + Al*Cos (pil80*15 . 0* (LST-tl) ) + TPAR 83 

& A2*Cos (pil80*30 . 0* (LST-t2) ) ) TPAR 84 


The SCALE parameter does not depend significantly on time of day, so the global 
average value is simply adjusted for latitude dependence by program step: 


C... Zonal average temperature scale height (km) vs. latitude 
SCALE = SCALE* (1.14 - 0 . 18*Cos (pil80*lat ) ) 


TPAR 87 
TPAR 88 


Following computation of the local (latitude and time-of-day-dependent) 
thermospheric parameters by subroutine Thermpar, the values are adjusted for long-term and 
short-term variability via the "ES" array of the Stewart model for the effects of seasonal 
variations in surface pressure ("DR"), the effects of dust storms ("DUST"), and the user- 
selected adjustment parameters ("deltaZF" for the height of the base of the thermosphere, 
"deltaTEX" for exospheric temperature, and "deltaTF" for temperature at the base of the 
thermosphere). The adjustments are accomplished in the "STEWART2" subroutine as 
follows: 
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ZF 


STW2 38 


= ZFO * EXP (ES ( 8 ) + ES(9)) + DR + DUST + deltaZF 
TINF = TINFO * EXP ( ES ( 2 ) + ES(3)) + deltaTEX STW2 46 

TF = TFO * EXP (ES { 8 ) + ES(9)) + deltaTF STW2 47 


The two ways to adjust thermospheric values are: (1) through input of “STDL”, which 
controls the long-term variability by the Stewart ES parameters [ES(2) for TINF and ES(8) 
for ZF and TF], see section 6, and (2) through input of adjustment parameters deltaZF, 
deltaTEX, and deltaTF, see section 6. Short-term variability by the Stewart ES parameters 
[ES(3) for TINF and ES(9) for ZF and TF] are set automatically in the program, during 
computation of the high and low density perturbation magnitudes. 

2.4 Revised Thermospheric Wind Calculations 

Wind components in the Mars-GRAM thermospheric height region are computed 
as follows: 


VISCFAC = 0 . 04*VISC/ ( 1 . 0E6*DENS*VLL**2 ) DSTP134a 

DENOM = CORIOL**2 + VISCFAC**2 DSTP158 

C... Viscous-corrected areostrophic wind components DSTP159 

EWWIND = (CORIOL*FUG - VISCFAC*FVG) /DENOM DSTP160 

NSWIND = (CORIOL*FVG + VISCFAC*FUG) /DENOM DSTP161 


where "VISC" is the coefficient of molecular viscosity, "DENS" is the atmospheric 
density, "VLL" is a vertical scale parameter, and "CORIOL" is the coriolis factor. At 
lower altitudes, where the viscosity effect is small, the eastward (northward) wind 
component ("EWWIND" or “NSWIND”) becomes equal to the eastward (northward) 
areostrophic wind (on Earth called the geostrophic wind) component ("FUG" or “FVG”). 

The areostrophic wind blows parallel to isobars (lines of constant pressure). At higher 
altitudes, where the viscosity effect becomes dominant, the winds take on a significant 
cross-isobar component, with the viscous-corrected wind blowing perpendicular to the 
isobars (toward the low pressure side). In the thermosphere, isotherms (lines of constant 
temperature) tend to align somewhat with the isobars. Therefore, the viscous wind model 
above has significant components along the isotherms at lower thermospheric altitudes and 
significant cross-iostherm components at higher altitudes. 

The only change in the viscous wind model for Mars-GRAM version 3.4 was to reduce 
the coefficient (line DSTP134a) from 1 to 0.04. In the original Stewart model, the pressure 
had extremely weak horizontal gradients (yielding weak areostrophic wind component 
estimates). With the realistic latitude-longitude variability for pressure that is now 
incorporated into the Mars-GRAM 3.4 thermospheric model, pressure gradients (and 
areostrophic winds) are much larger in the thermosphere and require the smaller coefficient on 
the viscosity term. The value of 0.04 was derived to reproduce as closely as possible the wind 
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data in the four MTGCM cases used to derive the parameterizations for the other 
thermospheric variables. 


8 



3. Results from the New Thermospheric Model 

Figures 3. 1 and 3.2 compare plots of exospheric temperature versus latitude and time 
of day from the Mars Thermospheric Circulation Global Model (MTGCM) case MGS98L (see 
table 2. 1 for characteristics of the model cases) with similar estimates from the new Mars- 
GRAM thermospheric model. Figures 3.3 and 3.4 show analogous comparisons between 
MTGCM case MGS97L and Mars-GRAM thermospheric temperatures. 

Figure 3.4 also illustrates the wind vectors from the new thermospheric model and 
shows the expected strong cross-isotherm components, discussed in Section 2.4. Strong 
cross-isotherm components are also a feature of the upper levels of the MTGCM output (e.g. 
Bougher et al., 1990, figure 7a; Bougher et al., 1993, figure 4a). 

Figures 3.5 and 3.6 illustrate the latitude versus time of day dependence of 
atmospheric density at a height intermediate between the base of the thermosphere and the 
exosphere (180 km), for the MTGCM case MGS98L. Figures 3.7 and 3.8 are comparable 
plots for MTGCM case MGS97L and Mars-GRAM new thermospheric model. 

Figure 3.9 shows the MTGCM case MGS98L values of temperature (TF) at the base of 
the thermosphere versus latitude and time of day. Figure 3.10 is a plot from Mars-GRAM for 
temperature at a height of 130 km (near the base of the thermosphere) and for the time 
corresponding to the MGS97L case. Figure 3.10 also shows the wind components derived 
from Mars-GRAM, and illustrates the expected significant along-iostherm components. 

Figure 3.1 1 illustrates the latitude and time of day dependence for the height of the base of the 
thermosphere (ZF) for MTGCM case MGS97L. 

Characteristics of the new thermospheric model in Mars-GRAM are to have a 
predominant wave-one (one wavelength per 360 degrees of longitude) component for 
exospheric temperature, but a predominant wave-two component, especially at low latitudes, 
for the temperature and height of the base of the thermosphere. These features can be seen as 
general characteristics of the MTGCM data (in figures 3.1, 3.3, 3.5, 3.9 and 3.1 1). Although 
Mars-GRAM cannot reproduce the effects of components with wave numbers higher than 
two, the effects are evident in the MTGCM results. General features of latitude versus time of 
day variation of the MTGCM results are considered to be adequately represented by the new 
Mars-GRAM thermospheric model. 

Figures 3.12 and 3.13 compare the height versus time of day behavior for MTGCM 
case MGS97L and Mars-GRAM. The general structure of the MTGCM results is reproduced 
in the Mars-GRAM plot. 
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4. The New Sub-Solar Longitude Model 


Regular publication of latitude and longitude of the sub-solar point on the Mars surface 
began in 1984. Consequently, when Mars-GRAM was originally developed in 1989, the time 
series of longitude for the sub-solar point on Mars was limited to the time period 1984-1988. 
Mars ephemeris data are now available in computer-readable form from the Jet Propulsion 
Laboratory (JPL). The available data include not only prior time periods, but also projections 
for a number of years into the future. A set on Mars longitudes for the sub-solar point for 
1984 through 2003 was constructed by making use of this JPL data (sent to us via electronic 
file transfer by Myles Standish of JPL). 

From the limited (1984 to 1988) time series, it was apparent that the sub-solar 
longitude has a significant contribution from periods longer than the Mars orbit period (687 
days). An empirical Fourier harmonic fit (with three harmonic terms) and a basic period of 
696 days was used in the original Mars-GRAM version. This procedure is retained through 
Mars-GRAM version 3.34. However, it was recently noticed (Penny Niles, Lockheed Martin) 
that for the years near the end of this century, the sub-solar longitude values from Mars- 
GRAM could be in error by several degrees. 

With the longer time series (1984 to 2003) a new harmonic fit showed that the sub- 
solar longitude variations are much better reproduced by a (three harmonic) Fourier series 
with two basic periods. The two periods in the new parameterization are 687 days (the 
fundamental period of the Mars orbit) and 777 days. The new relations reproduce sub-solar 
longitude values within a few tenths of a degree over the data series. This accuracy is 
comparable to frequently -used methods to compute the solar position from Earth using the day 
of the year as input. For leap year versus non-leap year variations, methods using day of the 
year input are accurate to a few tenths of a degree. 

Details of the new calculations for sub-solar longitude are seen in program changes for 
the ORBIT subroutine in Appendix A. 
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5. New Climate Adjustment Parameters 

Mars-GRAM was developed from parameterizations to atmospheric data observed by 
the Mariner (orbiter) and Viking (orbiter and lander) missions and represented the Mars 
atmosphere during the 1970's. Recent observations by Clancy et al. (1990) indicate, in 
response to atmospheric cooling as the dusty atmosphere has cleared in the last two decades, 
that current Mars temperature profiles are distinctly cooler than those observed in the Viking 
era. Comparisons between Mars-GRAM 3.34 mid-latitude average temperature profiles with 
recent data from Clancy (provided by Rich Zurek, private communication), indicate about 20 
K cooling in the 40-50 km height range, about 15 K at 20-30 km, but little change in the 5-10 
km altitude region. This temperature change could have significant effects on atmospheric 
density at high altitudes, and, as pointed out by Clancy et al. (1990), is quite important in 
planning for Mars missions that involve aerobraking. 

The temperature profiles in Mars-GRAM are built from parameterizations that 
estimate temperatures at significant levels (0, 5, 15, 30, 50 and 75-km altitudes). To allow 
user-controlled input incorporating the effects of temperature changes since the Viking era, a 
set of climate factor values may be input to Mars-GRAM version 3.4. These factors, specified 
by values of the climate adjustment factors (CF array), are multiplicative factors that alter 
temperatures derived from the original parameterizations in Mars-GRAM. For example, if all 
temperature values (from surface to 75 km) are to be increased by 10%, use values of 1 . 1 for 
all CF array elements. For a decrease of all temperatures by 10% use CF values of 0.9. 
Separate CF values may be selected at each significant level altitude (see section 6). 

After the effects of the climate factor array are included in the temperature estimates, 
effects on the pressure and density variations with altitude are then computed by applying the 
hydrostatic balance equation and perfect gas law to the temperature profile values. 

Together with adjustment parameters for the thermospheric variables (deltaTEX, 
deltaTF, and deltaZF), the climate adjustment factors now allow maximum flexibility in 
adjusting Mars-GRAM results to reproduce the temperature, pressure, and density expected 
under any specific simulation case. Appendix A lists all changes in the program code required 
to incorporate the new climate change factors. 
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6. Revised Program Execution Procedures 


Operation of the interactive form of Mars-GRAM is unchanged. Note however, that 
input values for the climate adjustment factors (CF array) and thermospheric adjustment 
parameters are not accepted by the interactive form (only default values of CF = 1 and 
deltaTEX = deltaTF = deltaZF = 0 are allowed in interactive mode). 

In the batch form of Mars-GRAM, input is provided via the NAMELIST file INPUT. 
The complete set of values for the reference INPUT file is as follows: 


$ INPUT 


LSTFL 

= 

•LIST’ , 

OUTFL 

= 

1 OUTPUT ' 

MONTH 

= 

7, 

MDAY 

SB 

20, 

MYEAR 

= 

76, 

NPOS 

= 

11, 

I HR 

= 

12, 

IMIN 

= 

30, 

SEC 

= 

0.0, 

ALSO 

= 

0.0, 

INTENS 

= 

0.0, 

RADMAX 

= 

0.0, 

DUSTLAT 

= 

0.0, 

DUSTLON 

= 

0.0, 

F107 


68.0, 

STDL 

= 

0.0, 

MODPERT 

= 

3, 

NR1 

= 

1001, 

NVARX 

= 

1, 

NVARY 

= 

0, 

LOG SCALE 

= 

0, 

FLAT 

= 

22.0, 

FLON 

= 

48.0, 

FHGT 

= 

-0.5, 

DELHGT 

= 

10.0, 

DELLAT 

= 

0.0, 

DELLON 

= 

0.0, 

DELTIME 

= 

0.0, 

CFO 

= 

1.0 

CF5 

= 

1.0 

CF15 

= 

1.0 

CF30 

= 

1.0 

CF50 

= 

1.0 

CF75 

= 

1.0 

deltaZF 

as 

0.0 

deltaTF 

= 

0.0 

deltaTEX 

= 

0.0 


$END 


List file name (CON for console listing) 

Output file name 
month of year 
day of month 

year (4-digit; 1970-2069 can be 2-digit) 

max # positions to evaluate (0 = read from TRAJDATA) 

GMT hour of day 

minute of hour 

second of minute (for initial position) 

starting Ls value (degrees) for dust storm (0 = none) 

dust storm intensity (0.0 - 3.0) 

max. radius (km) of dust storm (0 or >10000 = global) 

latitude (deg) for center of dust storm 

West longitude (deg) of center of dust storm 

10.7 cm solar flux (10**-22 W/cm**2 at 1 AU) 

std. dev. for thermosphere variation (-3.0 to +3.0) 

perturbation model; l=random, 2=wave, 3=both 

starting random number (0 < NR1 < 30000) 

x-code for plotable output (l=hgt above ref. ellipse) 

y-code for 2-D plotable output (0 for 1-D plots) 

0=linear scale, l=log scale, 2=COSPAR % deviations 

initial latitude (N positive) , degrees 

initial longitude (West positive) , degrees 

initial height (km), above ref. ellipse 

height increment (km) between steps 

latitude increment (deg) between steps 

West longitude increment (deg) between steps 

time increment (sec) between steps 

climate adjustment factor at surface 

climate adjustment factor at 5 km 

climate adjustment factor at 15 km 

climate adjustment factor at 30 km 

climate adjustment factor at 50 km 

climate adjustment factor at 75 km 

adjustment for base of thermosphere (km) 

adjustment for temperature at height ZF (K) 

adjustment for exospheric temperature (K) 


The newly required inputs are CFO through CF75 (climate adjustment factors at significant 
levels 0, 5, 15, 30, 50 and 75 km) and the thermospheric adjustment parameters (deltaZF for 
base height of the thermosphere, deltaTF for temperature of the base of the thermosphere, and 
deltaTEX for the exospheric temperature). 

A new feature of Mars-GRAM is that values of the input parameters are needed as 
input into the NAMELIST INPUT files only if they differ from the default values listed above. 
Thus, the output for the above reference case is generated by using an empty INPUT file, 
namely, 
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$ INPUT 
$END 


The reference LIST output file for the reference input case above is given in Appendix 
B and contains values for latitude and longitude of the sub-solar point, the heliocentric 
distance (Mars orbital radius), and if the altitude is above 75 km, the local exospheric 
temperature and the local height and temperature of the base of the thermosphere. 
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10.7 cm Solar Flux at Mars 

r 

Figure 2.1. Global mean exospheric temperature (K), versus 10.7 cm solar flux. Data are from Table 2 1 
The regression line is equation (2.1). < 
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Heliocentric Distance to Mars, AU 

Figure 2.2. Global mean height (km) of the base of the thermosphere, versus heliocentric distance. Data 
are from Table 2.1. The regression line is equation (2.2). 
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Figure 2.3. Global mean temperature (K) at the base of the thermosphere, versus 10.7 cm solar flux. 
Data are from Table 2.1. The regression line is equation (2.3). 
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Figure 2.4. Global mean vertical thermospheric temperature scale (km), versus 10.7 cm solar flux. 
Data are from Table 2.1. The regression line is equation (2.4). 



MTGCM Case MGS98L. 2 
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Local Solar Time, hours past noon 

Figure 3.1. MTGCM model exospheric temperature (K) versus latitude and time of day for model 
MGS98L. See Table 2.1 for values of model case parameters. 
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Figure 3.2. Mars-GRAM model exospheric temperature (K) versus latitude and time of day for model 
case MGS98L. See Table 2.1 for values of model case parameters. 
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Local Solar Time, hours past noon 

Figure 3.3. MTGCM model exospheric temperature (K) versus latitude and time of day for model 
MGS97L. See Table 2.1 for values of model case parameters. 
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Local Solar Time, Mars hours 

Figure 3.4. Mars-GRAM model exospheric temperature (K) versus latitude and time of day for model 
case MGS97L. See Table 2.1 for values of model case parameters. Wind speed vectors 
range from 0 to 129 m/s. 



MTGCM MGS98L Case, 180 km Height 
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Local Solar Time, Mars hours 

Figure 3.5. MTGCM model density (log-base-10, kg/m 3 ) at 180 km height, versus latitude and time of 
day for model case MGS98L. See Table 2.1 for values of model case parameters. 




Mars-GRAM MGS98L Case, 180 km Height 
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Local Solar Time, Mars hours 

Figure 3.6. Mars-GRAM model density (log-base-10, kg/m^) at 180 km height, versus latitude and 
of day for model case MGS98L. See Table 2.1 for values of model case parameters. 
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See Table 2.1 for values of model case parameters. 




MarsGRAM MGS97L Case, 180 km Height 



26 


Local Solar Time, Mars hours 

Figure 3.8. Mars-GRAM model density (log-base- 10, kg/m 3 ) at 180 km height, versus latitude and time of 
day for model case MGS97L. See Table 2.1 for values of model case parameters. 
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Local Solar Time, Mars hours 

Figure 3.9. MTGCM model temperature (IQ at base of the thermosphere, versus latitude and time of day 
for model case MGS98L. See Table 2.1 for values of model case parameters. 
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Local Solar Time, Mars hours 

Figure 3.10. Mars-GRAM model temperature (K) at 130 km height, versus latitude and time of day for 
model case MGS97L. See Table 2.1 for values of model case parameters. Wind speed vectors 
range from 0 to 188 m/s. 



MTGCM Case MGS97L, ZF (km) 
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Local Solar Time, Mars hours 

Figure 3.11. MTGCM model height (km) of the base of the thermosphere, versus latitude and time of day 
for model case MGS97L. See Table 2.1 for values of model case parameters. 






MTGCM Case MGS97L, Latitude 62.5 N 
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Figure 3.12. MTGCM model temperature (K) , versus height and time of day at latitude 62.5 degre 
for model case MGS97L. See Table 2.1 for values of model case parameters. 
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Local Solar Time, Mars hours 

Figure 3.13. Mars-GRAM model temperature (K) , versus height and time of day at latitude 62.5 degrees North, 
for model case MGS97L. See Table 2.1 for values of model case parameters. 




Appendix A - Updates in the Mars-GRAM Program Code 


1. Comparing MARSGRAM.FOR for Mars-GRAM 3.34 and Mars-GRAM 3.4 


The official version number is changed. 


***** Mars-GRAM 3 . 34 \MARSGRAM. FOR 

C... Program Mars-GRAM Interactive version 3.34 - November 1, 1995 
C 

***** Mars-GRAM 3 . 4 \MARSGRAM . FOR 

C... Program Mars-GRAM Interactive version 3.4 - April 1, 1996 
C 

★ * * ★ * 


MARS 1 
MARS 2 

MARS 1 
MARS 2 


Output of the official version number to the LIST file is changed. 


***** Mars-GRAM 3 . 3 4 \MARSGRAM . FOR 

C MARS 10 

1 Format ( ' Mars-GRAM Interactive version 3.34 - November 1, 1995') MARS 11 

C MARS 12 

***** Mars-GRAM 3 . 4 \MARSGRAM . FOR 

C MARS 10 

1 Format (' Mars-GRAM Interactive version 3.4 - April 1, 1996') MARS 11 

C MARS 12 

* * * * ★ 


Climate factors (CF) and thermospheric adjustment parameters (deltaZF, deltaTF, and 
deltaTEX) are added to the common DATACOM. 


***** Mars-GRAM 3.34 \MARSGRAM . FOR 

Sc logscale, dustlat , dustlon, dusthgt, radmax, Rref ,modpert, alsO , 
Sc intens , iuO , iup, maxf iles 
COMMON /FILENAME/lstfl, outfl 
***** Mars-GRAM 3 . 4 \MARSGRAM . FOR 

Sc logscale, dustlat , dustlon, dusthgt , radmax, Rref , modpert , alsO , 
Sc intens , iuO , iup , maxf iles , CF ( 0 : 5 ) , deltaZF, deltaTF, deltaTEX 
COMMON /FILENAME/ Is tfl, outfl 


MARS 22 
MARS 23 
MARS 24 

MARS 22 
MARS 23 
MARS 24 


The climate factors (CF) are initialized to their default values of 1.0. The thermospheric 
adjustment parameters (deltaZF, deltaTF, and deltaTEX) are initialized to their default values of 0.0. 
User selected values of these parameters can be used in the batch form of Mars-GRAM. In the 
interactive form of Mars-GRAM, these parameters cannot be changed from their default values. 


***** Mars-GRAM 3.34 \MARSGRAM . FOR 



Write ( iuO , 10 ) 





MARS 

44 

10 

format { 1 Enter name for LIST 

file 

(CON 

for console 

listing) : ' ) 

MARS 

45 

***** 

Mars-GRAM 3 . 4 \MARSGRAM . FOR 








Write ( iuO , 10 ) 





MARS 

44 

C. . . 

Set CF climate factors to 1.0 

{ .ne 

. 1 available in 

Batch version) 

MARS 

44a 


Do 5 i = 0 , 5 





MARS 

44b 


CF ( i ) = 1.0 





MARS 

44c 

5 

Continue 





MARS 

44d 

C. . . 

Set deltaZF, deltaTF & deltaTEX to 

o 

o 

(available in Batch only) 

MARS 

44e 


deltaZF =0.0 





MARS 

44f 


deltaTF = 0.0 





MARS 

44g 


deltaTEX =0.0 





MARS 

44h 

10 

format ( ' Enter name for LIST 

file 

(CON 

for console 

listing) : ’ ) 

MARS 

45 


***** 
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2. Comparing MARSGRMB.FOR for Mars-GRAM 3.34 and Mars-GRAM 3.4 


The official version number is changed. 

***** Mars-GRAM 3.34 \MARSGRMB . FOR 

C... Program Mars-GRAM Batch version 3.34 - November 1, 1995 
C 

***** Mars-GRAM 3 . 4XMARSGRMB .FOR 

C. . . Program Mars-GRAM Batch version 3.4 - April 1, 1996 
C 

***** 


MARB 1 
MARB 2 

MARB 1 
MARB 2 


The default value of solar flux at 1AU (F107) is changed to 68, from 185. The 1AU value of 
F107 on the day of the Viking 1 landing (July 20, 1976) was 68. 

***** Mars-GRAM 3.34 \MARSGRMB . FOR 

C DUSTLON =0.0, ! West longitude (deg) of center of dust storm 

C F107 = 185.0, I 10.7 cm solar flux (10**-22 W/cm**2 at 1 AU) 

C STDL =0.0, ! std. dev. for thermosphere variation (-3.0 

***** Mars-GRAM 3 . 4 \MARSGRMB . FOR 

C DUSTLON =0.0, ! West longitude (deg) of center of dust storm 

C F107 =68.0, ! 10.7 cm solar flux (10**-22 W/cm**2 at 1 AU) 

C STDL =0.0, ! std. dev. for thermosphere variation (-3.0 

***** 


MARB 35 
MARB 36 
MARB 37 

MARB 35 
MARB 36 
MARB 37 


The comment about the COSPAR model comparison (LOGSCALE = 2) is added. 


***** Mars-GRAM 3.34 \MARSGRMB . FOR 
C ! 

C LOGSCALE =0, ! 1 for log-base- 10 scale plots, 

C ! 

C FLAT =22.0, ! initial latitude (N positive), 

***** Mars-GRAM 3 . 4\MARSGRMB . FOR 


plots) 
0 for linear 
scale 

degrees 


C 

C LOGSCALE = 0, 

C 

C FLAT = 22.0, 

★ * ★ ★ ★ 


! plots) 

! 0=regular linear scale, l=log-base-10 scale, 

! 2=percentage deviations from COSPAR model 
! initial latitude (N positive) , degrees 


MARB 44 
MARB 45 
MARB 46 
MARB 47 

MARB 44 
MARB 45 
MARB 46 
MARB 47 


Descriptions are added to the NAMELIST input for the new climate adjustment factors (CF) 
and thermospheric adjustment parameters (deltaZF, deltaTF, and deltaTEX). 


***** Mars-GRAM 3.34 \MARSGRMB . FOR 


c 

DELTXME 

St 

o 

o 

! time increment (sec) between 

steps 

MARB 

53 

c 

$END 





MARB 

54 

***** Mars-GRAM 3. 

4 \MARSGRMB . FOR 




c 

DELTIME 

= 

0.0, 

! time increment (sec) between 

steps 

MARB 

53 

c 

CFO 

= 

1.0 

! climate adjustment factor at 

surface 

MARB 

53a 

c 

CF5 

= 

1.0 

! climate adjustment factor at 

5 km 

MARB 

53b 

c 

CF15 

= 

1.0 

! climate adjustment factor at 

15 km 

MARB 

53c 

c 

CF30 

- 

1.0 

! climate adjustment factor at 

30 km 

MARB 

53d 

c 

CF50 

= 

1.0 

! climate adjustment factor at 

50 km 

MARB 

53e 

c 

CF7 5 

= 

1.0 

! climate adjustment factor at 

75 km 

MARB 

53 f 

c 

deltaZF 

= 

0.0 

! adjustment for base of thermosphere (km) 

MARB 

53g 

c 

deltaTF 

= 

0.0 

! adjustment for temperature at height ZF (K) 

MARB 

53h 

c 

deltaTEX 

= 

0.0 

! adjustment for exospheric temperature (K) 

MARB 

53 i 

c 






MARB 

53j 

c 

$END 





MARB 

54 


★ ★ ★* ★ 
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3. Comparing MARSSUBS.FOR for Mars-GRAM 3.34 and Mars-GRAM 3.4 


The official version number is changed. 


***** Mars-GRAM 3.34 \MARS SUBS . FOR 

C Subroutines and Functions for Interactive and Batch versions of ALBL 1 

C Mars-GRAM version 3.34 - November 1, 1995 ALBL 2 

c ALBL 3 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

C Subroutines and Functions for Interactive and Batch versions of ALBL 1 

C Mars-GRAM version 3.4 - April 1, 1996 ALBL 2 

C ALBL 3 

* ★ ★ * * 


The climate adjustment factors (CF) and the thermospheric adjustment parameters 
(deltaZF, deltaTF, and deltaTEX) are added to the argument list for subroutine ATMOS2. The local 
exospheric temperature (Texos) and local temperature of the base of the thermosphere (Tbase) are 
added as output arguments, to pass these values for output to the LIST file. 

***** Mars-GRAM 3.34 \MARSSUBS . FOR 

& dus tM , dus t A , H , TEMP , DENST , UPFCTR , LWFCTR , PRES , RSC , ZO , TMAX , TMIN , 

Sc TAVG , Brunt f , densur f , t factor , ZF, alsO, intens, iuO) 

REAL LWFCTR , MARSAU , INTENS 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

Sc dus tM , dus tA , H , TEMP , DENST , UPFCTR , LWFCTR , PRES , RSC , ZO , TMAX , TMIN , 

Sc TAVG/ Brunt f , densur f , tf actor , ZF # alsO , intens # iuO , CF, deltaZF , 

Sc deltaTF, deltaTEX, Texos /Tbase) 

REAL LWFCTR , MARSAU , INTENS 

★ * ★ ★ * 


ATM2 2 
ATM2 3 
ATM2 4 

ATM2 2 
ATM2 3 
ATM2 3a 
ATM2 4 


Comments about the input/output argument list for ATMOS2 are updated, to reflect the 
new argument variables. 


***** Mars-GRAM 3 . 34XMARSSUBS . FOR 



COMMON 

/THERM/F107, stdl 

ATM2 

5 


Dimension ES(0:11) 

ATM2 

6 

c 



ATM2 

7 

c 

CHGT 

SPACECRAFT HEIGHT ABOVE REFERENCE SURFACE (KM) (INPUT) 

ATM2 

8 

C 

CLAT 

APACECRAFT LATITUDE (DEGREES) (INPUT) 

ATM2 

9 

C 

CLON 

WEST LONGITUDE OF SPACECRAFT (DEGREES) (INPUT) 

ATM2 

10 

C 

MARSAU 

MARS ORBITAL RADIUS (AU) (INPUT) 

ATM2 

11 

C 

SUNLAT 

AREOCENTRIC LATITUDE OF SUN (DEGREES) (INPUT) 

ATM2 

12 

c 

SUNLON 

MARS WEST LONGITUDE OF SUN (DEGREES) (INPUT) 

ATM2 

13 

C 

ALS 

AREOCENTRIC LONGITUDE OF SUN ORBIT (INPUT) 

ATM2 

14 

c 

DATE 

JULIAN DATE (INPUT) 

ATM2 

15 

c 

dustM 

dust storm magnitude for average T and p effect 

ATM2 

16 

c 

(1 = full magnitude, 0 = no dust storm) (INPUT) 

ATM2 

17 

C 

dustA 

dust storm magnitude for daily amplitude T and p effect 

ATM2 

18 

C 

(1 = full magnitude, 0 = no dust storm) (INPUT) 

ATM2 

19 

c 

H 

SCALE HEIGHT AT SPACECRAFT POSITION (KM) (OUTPUT) 

ATM2 

20 

c 

TEMP 

TEMPERATURE AT SPACECRAFT POSITION (K) (OUTPUT) 

ATM2 

21 

c 

DENST 

MASS DENSITY AT SPACECRAFT POSITION (KG/M**3) (OUTPUT) 

ATM2 

22 

c 

UPFCTR 

UPPER DEVIATION FACTOR ON MASS DENSITY (OUTPUT) 

ATM2 

23 

c 

LWFCTR 

LOWER DEVIATION FACTOR ON MASS DENSITY (OUTPUT) 

ATM2 

24 

c 

PRES 

PRESSURE AT SPACECRAFT POSITION (N/M**2) (OUTPUT) 

ATM2 

25 

c 

RSC 

AREOCENTRIC RADIUS TO SPACECRAFT (KM) (OUTPUT) 

ATM2 

26 

c 

ZO 

LOCAL TERRAIN HEIGHT RELATIVE TO REFERENCE ELLIPSOID 

ATM2 

27 

c 

(KM) (INPUT) 

ATM2 

28 

c 

TMAX 

Daily maximum surface temperature (K) (OUTPUT) 

ATM2 

29 

c 

TMIN 

Daily minimum surface temperature (K) (OUTPUT) 

ATM2 

30 

c 

TAVG 

Daily average surface temperature (K) (OUTPUT) 

ATM2 

31 

c 

Bruntf 

Brunt-Vaisala frequency = Sqrt ( (g/T) * (dT/dZ + g/Cp) ) 

ATM2 

32 

c 

(OUTPUT) 

ATM2 

33 
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C densurf Density at surface (kg/m**3) (OUTPUT) 

C tfactor Density perturbation factor at base of thermosphere 

C (OUTPUT) 

C 2F Height of base of thermosphere (OUTPUT) 

C alsO Initial Ls for dust storm 

C INTENS Dust storm intensity (0. 0-3.0) 

C iuO unit number for messages 

C 

***** Mars -GRAM 3 . 4XMARSSUBS .FOR 
COMMON /THERM/F107 , stdl 
Dimension ES ( 0 : 1 1 ) , CF ( 0 : 5 ) 

C 

C CHGT SPACECRAFT HEIGHT ABOVE REFERENCE SURFACE (KM) (INPUT) 

C CLAT SPACECRAFT LATITUDE (DEGREES) (INPUT) 

C CLON WEST LONGITUDE OF SPACECRAFT (DEGREES) (INPUT) 

C MARSAU MARS ORBITAL RADIUS (AU) (INPUT) 

C SUNLAT AREOCENTRIC LATITUDE OF SUN (DEGREES) (INPUT) 

C SUNLON MARS WEST LONGITUDE OF SUN (DEGREES) (INPUT) 

C ALS AREOCENTRIC LONGITUDE OF SUN ORBIT (INPUT) 

C DATE JULIAN DATE (INPUT) 

C dustM dust storm magnitude for average T and p effect 

C (1 = full magnitude, 0 = no dust storm) (INPUT) 

C dustA dust storm magnitude for daily amplitude T and p effect 

C (1 = full magnitude, 0 = no dust storm) (INPUT) 

C H SCALE HEIGHT AT SPACECRAFT POSITION (KM) (OUTPUT) 

C TEMP TEMPERATURE AT SPACECRAFT POSITION (K) (OUTPUT) 

C DENST MASS DENSITY AT SPACECRAFT POSITION (KG/M**3) (OUTPUT) 

C UPFCTR UPPER DEVIATION FACTOR ON MASS DENSITY (OUTPUT) 

C LWFCTR LOWER DEVIATION FACTOR ON MASS DENSITY (OUTPUT) 

C PRES PRESSURE AT SPACECRAFT POSITION (N/M**2) (OUTPUT) 

C RSC AREOCENTRIC RADIUS TO SPACECRAFT (KM) (OUTPUT) 

C 20 LOCAL TERRAIN HEIGHT RELATIVE TO REFERENCE ELLIPSOID 

C (KM) (INPUT) 

C TMAX Daily maximum surface temperature (K) (OUTPUT) 

C TMIN Daily minimum surface temperature (K) (OUTPUT) 

C TAVG Daily average surface temperature (K) (OUTPUT) 

C Bruntf Brunt-Vaisala frequency = Sqrt ( (g/T) * (dT/d2 + g/Cp) ) 

C (OUTPUT) 

C densurf Density at surface (kg/m**3) (OUTPUT) 

C tfactor Density perturbation factor at base of thermosphere 

C (OUTPUT) 

C 2F Local height of base of thermosphere (OUTPUT) 

C alsO Initial Ls for dust storm 

C INTENS Dust storm intensity (0. 0-3.0) 

C iuO unit number for messages 

C CF climate adjustment factor array for temperature levels 

C delta2F adjustment in height (km) of base of thermosphere 

C deltaTF adjustment in temperature (K) at height ZF 

C deltaTEX adjustment in exospheric temperature (K) 

C Texos local exospheric temperature (K) (OUTPUT) 

C Tbase local temperature for base of exosphere (K) (OUTPUT) 

C 

* ★ *r ★ ★ 


ATM 2 34 
ATM2 35 
ATM2 36 
ATM2 37 
ATM2 38 
ATM2 39 
ATM2 40 
ATM2 41 

ATM2 5 
ATM2 6 
ATM2 7 
ATM2 8 
ATM2 9 
ATM2 10 
ATM2 11 
ATM2 12 
ATM2 13 
ATM2 14 
ATM2 15 
ATM2 16 
ATM2 17 
ATM2 18 
ATM2 19 
ATM2 20 
ATM2 21 
ATM2 22 
ATM2 23 
ATM2 24 
ATM2 25 
ATM2 26 
ATM2 27 
ATM2 28 
ATM2 29 
ATM2 30 
ATM2 31 
ATM2 32 
ATM2 33 
ATM2 34 
ATM2 35 
ATM2 36 
ATM2 37 
ATM2 38 
ATM2 39 
ATM2 40 
ATM2 40a 
ATM2 40b 
ATM2 40c 
ATM2 40d 
ATM2 40e 
ATM2 40f 
ATM2 41 


The climate adjustment factors (CF) are added to the argument list of the subroutine 
TEMPS. 


***** Mars-GRAM 3.34 \MARSSUBS . FOR 

C. . . Evaluate temperatures at significant levels ATM2 68 

Call Temps (TObar, Tsurf, gam, CLAT, dlon, dustM, dustA, ALS, T) ATM2 69 

C... Evaluate surface pressure, N/m**2 ATM2 70 

***** Mars-GRAM 3 .4\MARSSUBS.FOR 

C. . . Evaluate temperatures at significant levels ATM2 68 

Call Temps (T0bar,Tsurf, gam, CLAT, dlon, dustM, dustA, ALS, T,CF) ATM2 69 

C. . . Evaluate surface pressure, N/m**2 ATM2 70 

★ ★ ★ ★ ★ 
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The local height of the base of the thermosphere (ZF) is computed, by calling the new 
subroutine Thermpar, and applying the Stewart (ES), pressure (DR), and dust corrections. The 
calling order for the Stewart2 subroutine is modified, to improve the program logic and run-time 
efficiency. This is the section that calls Stewart2 for use in the interpolation between 75 km altitude 
and the base of the thermosphere. The temperatures at the base of the thermosphere (Tbase) and in 
the exosphere (Texos) are passed for output to the LIST file. 


★ *★ 

c. . 

c. . 
c. . 


c. . 
c. . 

c. . 

★ ★ 

c. . . 

C. . . 

c. . . 


c. 


** Mars-GRAM 3 . 34\MARSSUBS . FOR 

Height of base of thermosphere 

ZF = 124.4 * <SMA / MARSAU) *EXP (ES { 8 ) ) + DR + DUST 
Shgt = height for call to thermosphere model 
Shgt = ZF if current height < ZF 
Shgt = CHGT 

If (CHGT .It. ZF) Shgt = ZF 

Evaluate atmospheric parameters at base of thermosphere 
Call Stewart 2 (MARSAU , ALS , CLAT , TIME , PRES , TEMP, DENST , 

& Shgt , Rs tar , H, AMF, 0.0, alsO , INTENS , iuO ) 

Evaluate density perturbation factor at base of thermosphere 
Call Stewart2 (MARSAU, ALS, CLAT, TIME, PRESHI , TEMPHI , 

Sc DENSHI , Shgt , Rs tar , HHI , AMHI ,1.0, als 0 , INTENS , iuO ) 

tfactor = DENSHI /DENST - 1.0 

Use stratosphere model if current height < ZF 
If (CHGT .It. ZF) then 
PF = PRES 

** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

Height of base of thermosphere 

Call Thermpar (MARSAU ,70., CLAT , TIME , SUNLAT , TINF0 , TF0 , ZF0 , SCALE ) 
ZF = ZF0*EXP (ES ( 8 ) ) + DR + DUST + deltaZF 
Use stratosphere model if current height < ZF 
If (CHGT .It. ZF ) then 
Shgt = ZF 

Evaluate thermospheric parameters at base of thermosphere 
Call S tewar 1 2 (MARSAU , ALS , CLAT , TIME , PRES , TEMP , DENST , Shgt , 
Rstar , H, AMF, 0.0, alsO , INTENS , iuO , sunlat, deltaZF, del taTF , 
del taTEX , TINF , TF ) 

Tbase = TF 
Texos = TINF 

Evaluate density perturbation factor at base of thermosphere 
Call Stewart2 (MARSAU, ALS, CLAT, TIME, PRESHI, TEMPHI, 

DENSHI , Shgt , Rstar , HHI , AMHI , 1 . 0 , alsO , INTENS , iuO , sunlat , 
deltaZF, deltaTF, deltaTEX, TINF, TF) 
tfactor * DENSHI / DENST - 1.0 
PF = PRES 


Sc 

Sc 


& 

Sc 


ATM2136 

ATM2137 

ATM2138 

ATM2139 

ATM2140 

ATM2141 

ATM2142 

ATM2143 

ATM2144 

ATM2145 

ATM2146 

ATM2147 

ATM2148 

ATM2149 

ATM2150 

ATM2151 

ATM2136 
ATM2 136a 
ATM2137 
ATM2149 
ATM2 149a 
ATM2149b 
ATM2149C 
ATM2 149d 
ATM2 1 4 9 e 
ATM2 1 4 9 f 
ATM2 1 4 9 g 
ATM2149h 
ATM2 1 4 9 i 
ATM2 1 4 9 j 
ATM2 149k 
ATM21491 
ATM2 149m 
ATM2151 


The calling order for the Stewart2 subroutine is modified, to improve the program logic and 
run-time efficiency. This is the section that calls Stewart2 within the thermosphere. The 
temperatures at the base of the thermosphere (Tbase) and in the exosphere (Texos) are passed for 
output to the LIST file. 


***** Mars-GRAM 3.34 \MARS SUBS . FOR 

Else ATM2163 

Call S tewar t2 (MARSAU, ALS, CLAT, TIME, PRESHI, ATM2164 

Sc TEMPHI , DENSHI , Shgt , Rs tar , HHI , AMHI ,1.0, alsO , INTENS , iuO ) ATM2 165 

Call S tewar t2 (MARSAU, ALS, CLAT, TIME, PRESLO, ATM2166 

Sc TEMPLO , DENSLO , Shgt , Rs tar , HLO , AMLO , - 1 . 0 , als 0 , INTENS , iuO ) ATM2 167 

UPFCTR = DENSHI /DENST ATM2168 

***** Mars-GRAM 3 . 4YMARSSUBS . FOR 

Else ATM2163 

Shgt = CHGT ATM2163a 

C... Evaluate thermospheric parameters at current height ATM2163b 

Call Stewart2 (MARSAU , ALS , CLAT , TIME , PRES , TEMP, DENST, ATM2163c 
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Sc 

& 


Sc 

& 

Sc 

Sc 


***** 


Shgt , Rs tar , H , AMF , 0 . 0 , als 0 , INTENS , iuO , sunlat , del taZF , del taTF , 
de 1 t aTEX , T INF , TF ) 

Texos = TINF 
Tbase = TF 

tfactor = DENSHI/DENST - 1.0 

Cal 1 Stewart2 (MARSAU , ALS , CLAT , TIME , PRESHI , TEMPHI , 

DENSHI , Shgt , Rs tar , HHI , AMHI ,1.0, alsO , INTENS , iuO , sunlat , 
deltaZF, del taTF , deltaTEX, TINF, TF) 

Call Stewart2 (MARSAU, ALS, CLAT, TIME, PRESLO, TEMPLO, 

DENSLO, Shgt, Rstar ,HLO, AMLO, -1.0, alsO, INTENS, iuO, sunlat, 
deltaZF, deltaTF , deltaTEX, TINF, TF) 

UPFCTR = DENSHI/DENST 


ATM2163d 

ATM2163e 

ATM2 1 6 3 f 

ATM2163g 

ATM2163h 

ATM2164 

ATM2165 

ATM2165a 

ATM2166 

ATM2167 

ATM2167a 

ATM2168 


Adjustment parameters CF, deltaZF, deltaTF, and deltaTEX are added to the common 
DATACOM. 


***** M ars-GRAM 3 . 34NMARSSUBS .FOR 

& logscale,dustlat,dustlon,dusthgt,radmax,Rref ,modpert, alsO , BLKD 

Sc intens , iuO , iup,maxf iles BLKD 

C... Unit numbers for messages (normally to screen) and to list output BLKD 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

S: logscale, dustlat , dustlon, dusthgt , radmax, Rref ,modpert, alsO , BLKD 

Sc intens , iuO , iup,maxf iles, CF (0 : 5) , deltaZF, deltaTF, deltaTEX BLKD 

C... Unit numbers for messages (normally to screen) and to list output BLKD 
***** 


8 

9 

10 

8 

9 

10 


Data statements are added, to establish default parameters. Note, the parameter maxfiles 
may need to be 16, for PC-based, 16-bit compilers. This will limit the number of output files that will 
be produced; maxfiles may be set to a larger value on UNIX systems and 32-bit PC compilers. 

***** Mars-GRAM 3.34 \MARSSUBS . FOR 
Data maxf iles/16/ 

C... Zurek wave perturbation parameter data values 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 
Data maxf iles/16/ 

C... Default = no dust storm 

Data AlsO , Intens , Radmax, Dustlat , Dustlon/5*0 . 0/ 

C. . . Default model parameters 

Data Modper t , NVARX , NVARY , LOGSCALE , NPOS /3, 1,0, 0,21/ 

C... Zurek wave perturbation parameter data values 
***** 


Adjustment parameters CF, deltaZF, deltaTF, and deltaTEX are added to the common 
DATACOM. 

***** Mars-GRAM 3 . 34XMARSSUBS . FOR 

& logscale, dustlat , dustlon, dusthgt , radmax, Rref, modper t, alsO, 

Sc intens, iuO , iup, maxf iles 
Real intens 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

Sc logscale, dustlat, dustlon, dusthgt , radmax, Rref , modper t, alsO , 

Sc intens, iuO , iup , maxf iles , CF (0:5) , deltaZF, deltaTF, deltaTEX 
Real intens 

***** 


DSTF 

11 

DSTF 

12 

DSTF 

13 

DSTF 

11 

DSTF 

12 

DSTF 

13 


BLKD 

18 

BLKD 

19 

BLKD 

18 

BLKD 

18a 

BLKD 

18b 

BLKD 

18c 

BLKD 

18d 

BLKD 

19 


Adjustment parameters CF, deltaZF, deltaTF, and deltaTEX are added to the common 
DATACOM. 


***** Mars-GRAM 3 . 34NMARSSUBS . FOR 

Sc logscale, dustlat , dustlon, dusthgt, radmax, Rref ,modpert, alsO , DSTP 9 

& intens , iuO , iup, maxf iles DSTP 10 
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COMMON /FILENAME/ lstfl.outfl 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

& logscale, dustlat , dustlon, dusthgt , radmax, Rref ,modpert, alsO, 
& intens, iuO , iup,maxf iles, CF (0 : 5) , deltaZF, deltaTF, deltaTEX 
COMMON /FILENAME/ lstfl.outfl 


DSTP 11 

DSTP 9 
DSTP 10 
DSTP 11 


The new input and output parameters (CF, deltaZF, deltaTF, deltaTEX, Texos, and Tbase) are 
added to the argument list in the call statement for subroutine ATMOS2. 


***** Mars-GRAM 3 . 34\MARSSUBS .FOR 

& dus tM , dus tA , HSCALE , TEMP , DENS , FACTHI , FACTLO , PRES , RSC , thgt , 

Sc TMAX, TMIN, TAVG , Brunt f , densurf , t factor , ZF, alsO, INTENS, iuO) 

C. . . CHGT = height above reference ellipsoid 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

& dus tM , dus tA , HSCALE , TEMP , DENS , FACTHI , FACTLO , PRES , RSC , thgt , 

Sc TMAX, TMIN, TAVG, Bruntf, densurf , tf actor, ZF, alsO, INTENS, iuO,CF, 
Sc deltaZF, deltaTF, deltaTEX, Texos , Tbase) 

Zbase = ZF 

C... CHGT = height above reference ellipsoid 


DSTP 41 
DSTP 42 
DSTP 43 

DSTP 41 
DSTP 42 
DSTP 42a 
DSTP 42b 
DSTP 43 


The new input and output parameters (CF, deltaZF, deltaTF, deltaTEX, Texos, and Tbase) are 
added to the argument list in the call statement for subroutine ATMOS2. 


***** Mars-GRAM 3.34 \MARSSUBS . FOR 

& dus tM , dus tA , HLONP , TLONP , DLONP , FHI , FLO , PLONP , RLONP , thgt , 

& TMAP , TMI P , TAVP , bvf , denO , t f ac , z t , als 0 , INTENS , iuO ) 

CLONM = CLON -2.5 
***** Mars-GRAM 3 . 4 \ MARS SUBS. FOR 

Sc dus tM , dus tA , HLONP , TLONP , DLONP , FHI , FLO , PLONP , RLONP , thgt , 

Sc TMAP , TMI P , TAVP , bvf , denO , tfac , zt , alsO , INTENS, iuO , CF, deltaZF, 
Sc deltaTF, deltaTEX, T INF, TF) 

CLONM = CLON - 2.5 


DSTP 87 
DSTP 88 
DSTP 89 

DSTP 87 
DSTP 88 
DSTP 88a 
DSTP 89 


The new input and output parameters (CF, deltaZF, deltaTF, deltaTEX, Texos, and Tbase) are 
added to the argument list in the call statement for subroutine ATMOS2. 


***** Mars-GRAM 3 . 34 \MARSSUBS . FOR 

Sc dus tM, dustA, HLONM, TLONM, DLONM, FHI , FLO , PLONM, RLONM, thgt , 

& TMAM, TMIM , TAVM, bvf , denO , tf ac , zt , alsO , INTENS, iuO) 

C... First two heights for area average height 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

Sc dus tM , dus tA , HLONM , TLONM , DLONM , FHI , FLO , PLONM , RLONM , thgt , 

Sc TMAM, TMIM, TAVM, bvf, denO, tfac, zt,als0, INTENS , iuO , CF, deltaZF, 
Sc deltaTF, deltaTEX, TINF,TF) 

C... First two heights for area average height 


DSTP 92 
DSTP 93 
DSTP 94 

DSTP 92 
DSTP 93 
DSTP 93a 
DSTP 94 


The new input and output parameters (CF, deltaZF, deltaTF, deltaTEX, Texos, and Tbase) are 
added to the argument list in the call statement for subroutine ATMOS2* 


***** Mars-GRAM 3 . 34NMARSSUBS .FOR 

Sc dustM, dus tA, HLATP, TLATP,DLATP, FHI, FLO, PLATP,RLATP, thgt, DSTP105 

Sc TMAP , TMIP , TAVP , bvf , denO , tfac , z t , alsO , INTENS , iuO ) DSTP10 6 

CLATM = CLAT - 2.5 DSTP1Q7 

***** Mars-GRAM 3 . 4\MARSSUBS .FOR 

& dustM, dustA, HLATP , TLATP , DLATP , FHI , FLO , PLATP , RLATP , thgt , DSTP 10 5 

Sc TMAP, TMIP, TAVP, bvf , denO, tfac, zt, alsO, INTENS, iuO,CF, deltaZF, DSTP106 

Sc deltaTF, deltaTEX, TINF,TF) DSTP106a 

CLATM = CLAT - 2.5 DSTP107 



The new input and output parameters (CF, deltaZF, deltaTF, deltaTEX, Texos, and Tbase) are 
added to the argument list in the call statement for subroutine ATMOS2. 


***** Mars-GRAM 3 . 34XMARSSUBS .FOR 

& dustM, dustA, HLATM, TLATM, DLATM, FHI , FLO, PLATM, RLATM, thgt , DSTP114 

& TMAM, TMIM, TAVM, bvf, denO , tf ac , zt , alsO , INTENS, iuO) DSTP115 

C... Area average height DSTP116 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

& dustM, dustA, HLATM, TLATM, DLATM, FHI, FLO, PLATM, RLATM, thgt, DSTP114 

& TMAM, TMIM, TAVM, bvf, denO , tfac , zt , alsO, INTENS, iuO,CF, deltaZF, DSTP115 

& deltaTF, deltaTEX, TINF,TF) DSTP115a 

C... Area average height DSTP116 

* -k ★ ★ ★ 


The coefficient on the viscosity term for the wind calculations is changed to 0.04 (from 1.0). 


***** Mars-GRAM 3 . 34XMARSSUBS .FOR 

VISC = BETA*TEMP**1.5/ (TEMP + SVAL) DSTP133 

VISCFAC = VISC/ (1.0E6*DENS*VLL**2) DSTP134 

C... FUG, FVG = Coriolis parameter times components of areostrophic DSTP135 
***** Mars-GRAM 3 . 4XMARSSUBS .FOR 

VISC = BETA*TEMP* *1.5/ ( TEMP + SVAL) DSTP133 

C. . . Use multiplier to scale wind magnitude DSTP134 

VISCFAC = 0 . 04*VISC/ ( 1 . 0E6*DENS*VLL**2 ) DSTP134a 

C. . . FUG, FVG = Coriolis parameter times components of areostrophic DSTP135 
★ * * * * 


The allowable upper and lower limits on the density perturbations are adjusted to better 
account for the stability limits that would be encountered. 


***** Mars-GRAM 3 . 34\MARSSUBS .FOR 

IF ( RHO . GE .0.0) DENSP = DENS*(1. + wave) + RHO*DPLUS DSTP237 

If (DENSP .It. 0.0) DENSP = 0.05*DENS DSTP238 

C. . . Standard deviation in random density perturbation (% of mean) , DSTP239 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

IF (RHO. GE. 0.0) DENSP = DENS*(1. + wave) + RHO*DPLUS DSTP237 

C. . . Re-check upper and lower bounds on density perturbations DSTP237a 

If (DENSP .It. 0.05*DENS) DENSP = 0.05*DENS DSTP238 

If (DENSP .gt. (l.+pertmax) *DENS) DENSP = ( 1 . +pertmax) *DENS DSTP238a 

C. . . Standard deviation in random density perturbation (% of mean) , DSTP239 

★ it * ★ ★ 


A re-check of the density perturbation values is done, by comparing with the stability limits. 


***** Mars-GRAM 3.34 \MARSSUBS . FOR 

SIGD = 50 . * (DENSHI-DENSLO) /DENS 
C. .. Add Zurek wave amplitudes to DENSHI, 
***** Mars-GRAM 3 . 4XMARSSUBS . FOR 

SIGD = 50 . * (DENSHI-DENSLO) /DENS 
C... Re-check SIGD, DENSHI and DENSLO 
If (SIGD . gt . 100 . *pertmax) Then 
SIGD = 100.*pertmax 
DENSHI = DENS* ( 1 . + pertmax) 

DENSLO = DENS/ ( 1 . + pertmax) 

Endif 

C... Add Zurek wave amplitudes to DENSHI, 
* ★ * 


DSTP241 

subtract from DENSLO DSTP242 

DSTP241 

DSTP241a 

DSTP241b 

DSTP241C 

DSTP241d 

DSTP241e 

DSTP241f 

subtract from DENSLO DSTP242 


The Mars orbital radius (MARSAU) and the solar longitude (SUNLON) are added to the 
output for the LIST file. 

***** Mars-GRAM 3 . 34 \MARSSUBS . FOR 

If (iup.gt.O ) Write (iup, 590) CSEC,CSEC/onesol,ALS,OHGT,OHGTS, DSTP248 
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& HSCALE, CLAT , CLON , SUNLON , TLOCAL DSTP249 

590 FORMAT ( * Time (rel. to TO) =’,F10.1,' sec. ( 1 , F8 . 3 , 1 sols) ' , DSTP250 

M ars -GRAM 3 . 4YMARSSUBS . FOR 

If (iup.gt .0) Write (iup, 590) CSEC, CSEC/onesol, ALS,OHGT, OHGTS, DSTP248 

& HSCALE , CLAT , CLON , SUNLAT , MARSAU , SUNLON , TLOCAL DSTP249 

590 FORMAT { ' Time (rel. to TO) =’,F10.1, 1 sec. ( ' , F8 . 3 , ’ sols) ' , DSTP250 

★ * ★ ★ ★ 


The formats are revised to accommodate the new output values. The line numbers are 
changed on several lines to allow for these changes. 


***** Mars -GRAM 3 . 34 \MARSSUBS . FOR 

Sc 'Scale Height = ' , F6.2,* km'/' Latitude = ' ,F8.3,' degrees ' , 
& ' West Longitude = *,F8.3,' degrees’/’ Sun Longitude = 

Sc F8 . 3 # ' deg. Local Time = ' ,F6.2, ' Mars hours') 

Call cospar ( OHGT f tcos , pcos , dcos ) 

If (dcos . le . 0 . 0) Then 
devlo = -99.9 
devav = -99.9 
devhi = -99.9 
Else 

devlo = 100 . * (DENSLO-dcos ) /dcos 
devav = 100 . * (DENS-dcos ) /dcos 
devhi = 100. * (DENSHI-dcos) /dcos 
End if 

I f ( iup . gt . 0 ) Wr i te ( iup ,600) TEMP , PRES , DENSLO , DENS , DENSHI , 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

Sc 'Scale Height = ',F6.2,' km'/' Latitude = 1 ,F8.3,' degrees ' , 

Sc ' West Longitude = ',F8.3,' degrees'/' Sun Latitude =',F9.2, 

Sc ' deg. Mars Orbital Radius =',F6.3,' AU'/' Sun Longitude =' 

Sc ,F8.2,' deg. Local Time = ',¥6.2,' Mars hours’) 

If ( iup . gt . 0 . and . OHGT . gt . 7 5 . 0 ) Write ( iup , 595 ) Texos , Tbase , Zbase 
595 Format (' Exospheric Temp. = ' ,F6.1,’ K',7X,'Tbase = ' ,F6.1, 

Sc ’ K ' , 3X, ' Zbase = \F6.1,’ km') 

Call cospar (OHGT, tcos , pcos , dcos ) 

If (dcos . le . 0 . 0) Then 
devlo = -99.9 
devav = -99.9 
devhi = -99.9 
Else 

devlo = 100 .* (DENSLO-dcos ) /dcos 
devav = 100 . * (DENS-dcos) /dcos 
devhi = 100 . * (DENSHI-dcos) /dcos 
End if 

I f ( iup . gt . 0 ) Wr i t e ( iup ,600) TEMP , PRES , DENSLO , DENS , DENSHI , 


DSTP252 

DSTP253 

DSTP254 

DSTP254b 

DSTP254C 

DSTP254d 

DSTP254e 

DSTP254 f 

DSTP254g 

DSTP2 54h 

DSTP2 54i 

DSTP2 54 j 

DSTP2 54k 

DSTP255 

DSTP252 

DSTP252a 

DSTP253 

DSTP254 

DSTP254a 

DSTP254b 

DSTP254c 

DSTP254d 

DSTP254e 

DSTP254f 

DSTP254g 

DSTP2 54h 

DSTP254i 

DSTP2 54 j 

DSTP2 54k 

DSTP2541 

DSTP2 54m 

DSTP255 


The output variable for NVARX=8 is changed to local solar time in hours (from hour angle 
in degrees). 


***** Mars-GRAM 3 . 34\MARSSUBS . FOR 
If (NVARX.EQ.7)VARX = ALS 
I f ( NVARX . EQ . 8 ) VARX = TLOCAL* 15. 
Alogdens = 0.0 

***** Mars-GRAM 3 . 4\MARSSUBS . FOR 
If (NVARX. EQ.7JVARX = ALS 
If (NVARX. EQ. 8) VARX = TLOCAL 
Alogdens = 0.0 


DSTP272 

DSTP273 

DSTP273a 

DSTP272 

DSTP273 

DSTP273a 


The output variable for NVARY=8 is changed to local solar time in hours (from hour angle 
in degrees). 
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u * 


***** 


Mars-GRAM 3 .34\MARSSUBS.FOR 
I f (NVARY . EQ . 7 ) VARY = ALS 
I f (NVARY . EQ . 8 ) VARY = TL0CAL*15. 
WRITE ( 2 1 , 7 9 5 ) VARX , VARY , SIGD 
***** Mars-GRAM 3 . 4XMARSSUBS .FOR 
If (NVARY. EQ. 7) VARY = ALS 
If (NVARY. EQ. 8) VARY = TLOCAL 
WRITE ( 2 1 , 7 9 5 ) VARX , VARY , SIGD 


DSTP324 

DSTP325 

DSTP326 

DSTP324 

DSTP325 

DSTP326 


The new coefficients are loaded for the 687 day plus 777 day sun longitude calculations 
(previously a single, 696 day, periodicity was assumed). 


***** Mars-GRAM 3 .34\MARSSUBS.FOR 

C... d0-d6 = Fourier coefficients for lonsun calculation 

data dO ,dl, d2 , d3 , d4, d5,d6/-3 . 5158255d0 , 6 . 5240771d0, -9.3284323d0, 

Sc 1.9052424d0, 2. 4887938d0, -0. 49182825 dO ,0.10 13 0513d- 1/ 

C... perl = 687 day period; per2 = 696 day period 
perl=6 . 87d2 /TWOPI 
per2=6 . 96d2/TWOPI 
PI180 = TWOPI/3 . 60d2 
***** Mars-GRAM 3 . 4XMARSSUBS .FOR 

C. . . d0-d6 = Fourier coefficients for lonsun calculation 

Data eO , el , e2 , e3 , e4 , e5 , e6/- . 61024174D-01 , 1 . 8014752D0 , . 35080142D0 , 
& . 1182 052 8D0 , - . 23869660D0 ,-.56879633D-01,-. 42787450D-01/ 

Data dO , dl , d2 , d3 , d4 , d5 , d6/-3 . 4766116DO, 10 . 199478D0 , 2 .7450878D0 , 

& -3 . 2726522D0 , - . 67435053D-01 , .52016400DO, .11903292D0/ 

C... perl = 687 day period; per2 = 777 day period 
perl=6.87d2 /TWOPI 
per2=7 . 77d2 /TWOPI 
PI180 = TWOPI/3. 60d2 


ORBT 12a 
ORBT 13 
ORBT 14 
ORBT 14a 
ORBT 15 
ORBT 16 
ORBT 17 

ORBT 12a 
ORBT 12b 
ORBT 12c 
ORBT 13 
ORBT 13a 
ORBT 14 
ORBT 15 
ORBT 16 
ORBT 17 


The co mm ent line is modified and the computation of the time variable TIME2 is moved. 


***** Mars-GRAM 3 .34XMARS SUBS. FOR 

xlon = lonO + 3.60d2*(date - 2 . 922d3) /perlon 
C. . . TIME1 , TIME2 = time variables for periods PERI and PER2 
TIME1=DATE/PER1 
TIME2=DATE/PER2 

... Lsubs calculation (REAL*8) 

**** Mars-GRAM 3 . 4\MARSSUBS .FOR 

xlon = lonO + 3.60d2*(date - 2 . 922d3 ) /perlon 
C . . . TIME1 = time variable for period PERI 
TIME1=DATE/PER1 

C... Lsubs calculation (REAL*8) 


ORBT 20 
ORBT 20a 
ORBT 21 
ORBT 22 
ORBT 22a 

ORBT 20 
ORBT 20a 
ORBT 21 
ORBT 22a 


The new, dual period, solar longitude model equations are evaluated. 


***** Mars-GRAM 3 . 34\MARSSUBS .FOR 
C... lonsun calculation (REAL*8) 

LON=DO+D1*DSIN (TIME2 ) +D2*DCOS (TIME2) +D3*DSIN (2 . 0d0*TIME2 ) 

Sc +d4*DCOS (2 . 0d0*time2 ) +d5*DSIN (3 . 0d0*time2 ) +d6*DCOS ( 3 . 0d0*time2 ) 
& +xlon 

C... Put Lsubs and lonsun into 0-360 degree range 
LS = DMOD (LS , 3 . 6D2 ) 

***** Mars-GRAM 3 . 4XMARSSUBS . FOR 
C... lonsun calculation (REAL* 8) 
time = (date-2 . 921d3 ) /perl 
time2 = (date-2 . 921d3 ) /per2 

xlon=D0+Dl*DSIN (TIME) +D2*DCOS (TIME) +D3*DSIN(2 . 0d0*TIME) 

Sc +d4*DCOS (2 . 0d0*time) +d5*DSIN (3 . 0d0*time) +d6*DCOS (3 . 0d0*time) 

Sc + xlon 

LON=E0+El*DSIN(time2) +E2*DCOS(time2) +E3*DSIN(2 .0d0*time2) 


ORBT 27a 
ORBT 28 
ORBT 29 
ORBT 30 
ORBT 30a 
ORBT 31 

ORBT 27a 
ORBT 27b 
ORBT 27c 
ORBT 28 
ORBT 29 
ORBT 30 
ORBT 30a 
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& +E4*DC0S (2 . 0d0*time2 ) +E5*DSIN (3 . OdO*time2 ) +E6*DC0S (3 . 0d0*time2 ) 
& + xlon 

C... Put Lsubs and lonsun into 0-360 degree range 
LS = DMOD(LS, 3 . 6D2 ) 

★ * * * ★ 


The new variables (sunlat, deltaZF, deltaTF, deltaTEX, TEMF, and TF) are added to the 
argument list for the call to subroutine Stewart2. 

***** Mars-GRAM 3 . 34 \MARSSUBS . FOR 

SUBROUTINE STEWART 2 ( RAUI , LSUN, LAT , LST, TOTAL PRZ , TZ, 

& TOTALMDZ , CHGT, RSTAR, H, MOLWTG, SIGMA, alsO, INTENS , iuO) 

C 

***** Mars-GRAM 3 . 4 \MARSSUBS - FOR 

SUBROUTINE STEWART 2 ( RAUI, LSUN, LAT, LST, TOTAL PRZ , TZ, 

& TOTALMDZ, CHGT, RSTAR, H, MOLWTG, SIGMA, alsO, INTENS, iuO , 

& sunlat, deltaZF, deltaTF, deltaTEX, T INF, TF) 

C 

★ ★ * * * 


A comment line is added to clarify that FBARR is for the Mars orbital position. 

***** Mars-GRAM 3.34 \MARSSUBS . FOR 
RAU = RAUI 

FBARR = FBAR / (RAU* *2) 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 
RAU = RAUI 

C... Convert solar 10.7 cm flux to Mars position 
FBARR = FBAR / (RAU**2) 

* * * ★ ★ 

The local values of the thermospheric parameters are evaluated by railing the new 
subroutine Thermpar. 

***** Mars-GRAM 3 . 34NMARSSUBS .FOR 
CALL PRSEAS (LSUN, LAT, PFAC) 

if (FLAG. gt.O) Write (iuO, *) ' RREF, RAU, GZ = ' , RREF, RAU, GZ 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

CALL PRSEAS (LSUN, LAT, PFAC) 

C... Evaluate the basic parameters for the thermosphere model 

Call Thermpar ( RAU , FBARR , LAT , LST , SUNLAT , TINFO , TFO , ZFO , SCALE ) 
if (FLAG. gt . 0) Write (iuO , * ) ' RREF, RAU, GZ = 1 , RREF , RAU , GZ 

* * ★ ★ * 

The modifications, for the Stewart adjustments (ES), pressure (DR), dust, and the new 
adjustment parameter deltaZF, are added to the value (ZFO) computed by Thermpar. 


***** Mars-GRAM 3 . 34\MARSSUBS . FOR 

DUST = DUST * EXP (ES ( 10 ) ) STW2 37 

ZF = (124.4 * SMA / RAU) * EXP(ES(8> + ES(9)) + DR + DUST STW2 38 

ZZF = CHGT — ZF STW2 39 

***** Mars-GRAM 3 . 4\MARSSUBS . FOR 

DUST = DUST * EXP (ES (10) ) STW2 37 

C... Height- of base of thermosphere STW2 37a 

ZF = ZFO * EXP(ES{8) + ES{9)) + DR + DUST + deltaZF STW2 38 

C... Height above base of thermosphere STW2 38a 

ZZF = CHGT — ZF STW2 3 9 


The modifications, for the Stewart adjustments (ES) and the new adjustment parameter 
deltaTEX, are added to the value (TINFO) computed by Thermpar. 


STW2 31 
STW2 32 

STW2 31 
STW2 31a 
STW2 31b 
STW2 32 


STW2 28 
STW2 29 

STW2 28 
STW2 28a 
STW2 29 


STW2 1 
STW2 2 
STW2 3 

STW2 1 
STW2 2 
STW2 2a 
STW2 3 


ORBT 30b 
ORBT 30c 
ORBT 30d 
ORBT 31 
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* o 


** **★ 


Mars-GRAM 3 . 34\MARSSUBS . FOR 



RF = RREF + ZF 

STW2 

45 


TINF = 4.11 * (11.0 + FBARR) * EXP(ES{2) + ES(3)) 

STW2 

46 


TF = (170.0 * SMA / RAU) * EXP(ES(8) + ES(9)) 

STW2 

47 

***** 

If (FLAG .GT. 0) then 
Mars-GRAM 3 . 4\MARSSUBS .FOR 

STW2 

48 


RF = RREF + ZF 

STW2 

45 

c. . . 

Exospheric temperature 

STW2 

45a 


TINF = TINF0 * EXP (ES (2 ) + ES(3)) + deltaTEX 

STW2 

46 

c. . . 

Temperature at base of thermosphere 

STW2 

46h 


TF = TF0 * EXP (ES ( 8 ) + ES(9)) + deltaTF 

STW2 

47 


If (FLAG .GT. 0) then 

STW2 

48 


★ ★Hr** 


The SCALE parameter is added as an argument for the THERMOS subroutine. 


***** Mars-GRAM 3 . 34XMARSSUBS .FOR 

CALL THERMOS (FLAG, ES, TINF, TF, LAT, LST, ZF,RF, ZZF, TOTAL PRZ , 
& TOTALNDZ, TZ, MOLWTG, PRZ, NDZ, MDZ , TOTALMDZ ,iuO) 

. . . SCALE HEIGHT, km 

**** Mars-GRAM 3 . 4NMARSSUBS . FOR 

CALL THERMOS (FLAG, ES, TINF, TF, LAT, LST, ZF , RF , ZZF, TOTALPRZ , 
& TOTALNDZ, TZ, MOLWTG, PRZ, NDZ, MDZ, TOTALMDZ ,iuO, SCALE) 

C... SCALE HEIGHT, km 


STW2 58 
STW2 59 
STW2 60 

STW2 58 
STW2 59 
STW2 60 


The SCALE parameter is added as an argument is the call statemnt for the THERMOS 
subroutine. 


***** Mars-GRAM 3 . 34XMARSSUBS .FOR 

SUBROUTINE THERMOS (FLAG, ES, TINF, TF, LAT, LST, ZF, RF, ZZF, 
Sc TOTALPRZ, TOTALNDZ, TZ, MOLWTG, PRZ, NDZ, MDZ, TOTALMDZ, iuO) 
INTEGER FLAG , iuO 

DIMENSION ES (0:11) , PRZ (0:11) ,NDZ(0:11) ,MDZ(0:11) 

***** Mars-GRAM 3 . 4\MARSSUBS . FOR 

SUBROUTINE THERMOS (FLAG, ES, TINF, TF, LAT, LST, ZF, RF, ZZF, 
Sc TOTALPRZ, TOTALNDZ, TZ, MOLWTG, PRZ, NDZ, MDZ, TOTALMDZ, iuO, 
Sc SCALE) 

INTEGER FLAG , iuO 

DIMENSION ES (0:11) ,PRZ(0:11) ,NDZ(0:11) ,MDZ(0:11) 


THRM 1 
THRM 2 
THRM 3 
THRM 4 

THRM 1 
THRM 2 
THRM 2a 
THRM 3 
THRM 4 


Computation of the SCALE parameter within the THERMOS subroutine is deleted, since its 
value now is input via the argument list 


***** Mars-GRAM 3 . 34\MARSSUBS . FOR 

C FO = PARAMETER IN EQUATION FOR ATOMIC OXYGEN CONTENT I THRM 39 

SCALE = TF / 9.20 THRM 40 

M(0) = 44.011 THRM 41 

***** Mars-GRAM 3 .4\MARSSUBS.FOR 

C. . . . FO = PARAMETER IN EQUATION FOR ATOMIC OXYGEN CONTENT I THRM 39 

M(0) = 44.011 THRM 41 

***** 


The climate adjustment factors (CF) are added to the argument list for the subroutine 
TEMPS. 


***** Mars-GRAM 3 . 34\MARSSUBS . FOR 

C THRM143 

Subroutine Temps (TObar, Tsurf, gam, alat, dlon,dustM, dustA, als , T) TMPS 1 

C TMPS 2 

***** Mars-GRAM 3 ,4\MARSSUBS.FOR 

C THRM143 


Subroutine Temps (TObar , Tsurf , gam, alat ,dlon, dustM, dustA, als , T, CF) TMPS 1 
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TMPS 2 


C 


The climate adjustment factors (CF) are added to the dimension statement for the 
subroutine TEMPS. 


***** Mars -GRAM 3 . 3 4 \MARSSUBS . FOR 

Dimension gam (5) , T (0 : 5) ,dz(5) ,z{5) , factor (5 ) , df actor (5 ) , af act { 5 ) TMPS 17 
Sc ,A25<0:6) , B2 5 ( 0 : 6 ) , C25 ( 0 : 6 ) TMPS 18 

C... Coefficients for T(25 km) correction to gammas TMPS 19 

***** Mars -GRAM 3 . 4 \MARSSUBS . FOR 

Dimension gam (5) , T (0 : 5 ) , dz (5 ) ,z(5) , factor ( 5 ) , df actor ( 5 ) , af act ( 5 ) TMPS 17 
& , A2 5(0:6) ,B2 5(0:6) ,C2 5(0:6) , CF (0:5) TMPS 18 

C... Coefficients for T{25 km) correction to gammas TMPS 19 

★ ★ * ★ * 


The temperature at the significant levels are evaluated, with the inclusion of a multiplicative 
climate adjustment factor (CF). 

***** Mars-GRAM 3 . 34 \MARSSUBS . FOR 
T ( 0 ) = Tsurf 

C. . . Re-evaluate lapse rates (deg. /km), based on actual temperatures 
***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 
T ( 0 ) = Tsurf 

C... Apply climate adjustment factors 
Do 45 i = 0,5 
T ( i ) = T(i) *CF<i) 

45 Continue 

C. . . Re-evaluate lapse rates (deg. /km), based on actual temperatures 
* * * * * 


TMPS138 

TMPS139 

TMPS138 
TMPS 13 8a 
TMPS 13 8b 
TMPS 13 8 c 
TMPS13 8d 
TMPS139 


The new subroutine Thermpar is added. Thermpar evaluates the local values of TINF0, the 
exospheric temperature (K), TFO, the temperature (K) at the base of the thermosphere, ZFO, the height (km) of 
the base of the thermosphere, and SCALE, the scale height (km) for the thermospheric temperature variations. 


***** Mars-GRAM 3 . 34\MARSSUBS . FOR 

c TMPS 144 

Subroutine Tsurface (sitla, sitlo , sunla, sunlo, als , au, Tsurf , TSRF 1 

***** Mars-GRAM 3 . 4 \MARSSUBS . FOR 

c TMPS144 

Subroutine Thermpar (RAU, FBARR, lat , LST, sunlat , TINF0 , TFO , ZFO , TPAR 1 

Sc SCALE ) TPAR 2 

REAL lat , LST , LATMAX TPAR 3 

Data LATMAX/ 2 5. 4/ TPAR 4 

C TPAR 5 

C TPAR 6 

C TPAR 7 

C. . . Thermospheric parameters, revised from the orginal Stewart TPAR 8 

C parameterizations : TPAR 9 

C SMA = 1.523691 TPAR 10 

C ZFO = 124.4 * (SMA/RAU) TPAR 11 

C TINF0 = 4.11 * (11.0 + FBARR) TPAR 12 

C TFO = 170.0 * (SMA/RAU) TPAR 13 

C SCALE = TFO / 9.20 TPAR 14 

C TPAR 15 

C The new parameterizations are based on four data sets from the TPAR 16 

C University of Arizona Mars Thermospheric Global Circulation TPAR 17 

C Model (MTGCM) , cases MGS97L , MGS98L, MANC00, and MGS 9 7E. For TPAR 18 

C a description of the MTGCM model and its output, see Bougher, TPAR 19 

C et al . , Journal of Geophysical Research, vol. 95 (B9) , pp . TPAR 20 

C 14,811 - 14,827, August 30, 1990. TPAR 21 

C TPAR 22 

C TPAR 23 

C TPAR 24 
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c. . . 


c. . . 


c. . . 
c. . . 


c 


Inputs : 
RAU 
FBARR 
lat 
LST 
sunlat 
Outputs : 
TINFO 
TFO 
ZFO 
SCALE 


orbital position radius (AU) 

10.7 cm solar flux at Mars position 
latitude for evaluation of parameters (degrees) 
local solar time (Mars hours) at evaluation point 
latitude of sun (degrees) 

Exospheric temperature (K) 

Temperature at base of thermosphere (K) 

Height of base of thermosphere (km) 

Scale height for temperature variations (km) 


TPAR 25 
TPAR 26 
TPAR 27 
TPAR 28 
TPAR 29 
TPAR 30 
TPAR 31 
TPAR 32 
TPAR 33 
TPAR 34 
TPAR 35 
TPAR 36 


Output values are un-corrected for Stewart (ES array) variations, 
pressure and dust effects. These factors are accounted for in 
the Stewart2 subroutine. Adjustment factors deltaTEX, deltaTF 
amd deltaZF are also added after computation of these values. 


Degrees to radians conversion factor 
pil80 = Atan (1 . ) /45 . 

Global mean exospheric temperature (K) versus 10.7 cm flux 
Tbar = 156.3 + 0.9427*FBARR 

Zonal average exospheric temperature (K) versus latitude 
Tavg = TbarMl. + 1 . 369E-4*sunlat*lat) 

Phase angles (hours) for local solar time variation 
tl = 13.2 - 0 . 00119*sunlat*lat 
t2 = 9.4 - 0 . 00231*sunlat*lat 

Amplitude factor for local solar time variation 
cphi = Cos (pil80* ( lat + sunlat) /(l. + LATMAX/90 . ) ) 

Exospheric temperature (K) versus local solar time 
TINFO = Tavg* ( 1 . + 0 . 22*cphi*Cos (pil80*15 . 0* (LST-tl) ) + 

& 0 . 04*cphi*Cos (pil80*30.0* (LST-t2) ) ) 

Global mean height of thermosphere base (km) 

Zbar = 197.94 - 49.058*RAU 

Latitude variation factor 

factlat = ( sunlat /LATMAX) * ( lat/ 77 . 5 ) **3 

Zonal average base height (km) versus latitude 

Zavg = Zbar + 4. 3* factlat 

Amplitudes for local solar time variation 
Al = 1.5 - Cos (pil80*4 . 0*lat) 

A2 = 2 . 3* (Cos (pil80* (lat + 0 . 5*sunlat) ) ) **3 
Phase angles (hours) for local solar time variation 
tl = 16.2 - (sunlat /LATMAX) * Atan (pi 180* 10 . 0*lat) 
t2 = 11.5 

Base height of thermosphere (km) versus local solar time 
ZFO = Zavg + Al*Cos (pil80* 15 . 0* (LST-tl ) ) + 

Sc A2*Cos (pil80*30.0* (LST-t2) ) 

Global mean temperature (K) at thermosphere base, versus FBARR 
Tbar = 113.7 + 0.5791*FBARR 

Zonal average temperature at thermosphere base (K) vs. latitude 
Tavg = Tbar*(l. + 0 . 186*factlat) 

Amplitudes for local solar time variation 
Al = 0.06 - 0 . 05*Cos (pil80*4 . 0*lat ) 

A2 = 0.1* (Cos (pil80* ( lat + 0 . 5*sunlat) ) ) **3 
Phase angles (hours) for local solar time variation 
tl = 17.5 - 2.5* (sunlat /LATMAX) *Atan (pil80*10. 0*lat ) 
t2 = 10.0 + 2.0*(lat/77.5)**2 

Thermosphere base temperature (K) versus local solar time 
TFO = Tavg* (1.0 + Al*Cos (pil80*15 . 0* (LST-tl) ) + 

& A2*Cos(pil80*30.0* (LST-t2) ) ) 

Global mean scale height (km) of thermospheric temperature 
SCALE = 8.38 + 0 . 09725*FBARR 

Zonal average temperature scale height (km) vs. latitude 
SCALE = SCALE* (1.14 - 0 . 18*Cos (pil80* lat ) ) 

Return 

End 


TPAR 37 
TPAR 38 
TPAR 39 
TPAR 40 
TPAR 41 
TPAR 42 
TPAR 43 
TPAR 44 
TPAR 45 
TPAR 46 
TPAR 47 
TPAR 48 
TPAR 49 
TPAR 50 
TPAR 51 
TPAR 52 
TPAR 53 
TPAR 54 
TPAR 55 
TPAR 56 
TPAR 57 
TPAR 58 
TPAR 59 
TPAR 60 
TPAR 61 
TPAR 62 
TPAR 63 
TPAR 64 
TPAR 65 
TPAR 66 
TPAR 67 
TPAR 68 
TPAR 69 
TPAR 70 
TPAR 71 
TPAR 72 
TPAR 73 
TPAR 74 
TPAR 75 
TPAR 76 
TPAR 77 
TPAR 78 
TPAR 79 
TPAR 
TPAR 
TPAR 
TPAR 
TPAR 
TPAR 
TPAR 86 
TPAR 87 
TPAR 88 
TPAR 89 
TPAR 90 
-TPAR 91 


80 

81 

82 

83 

84 

85 
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Subroutine Tsurface (sitla,sitlo, sunla , sunlo, als , au, Tsurf 


TSRF 1 


★ ★ * * * 


4. Comparing SETUP.FOR for Mars-GRAM 3.34 and Mars-GRAM 3.4 
The official version number is changed. 

***** Mars-GRAM 3 . 34\SETUP . FOR 


u SETU 3 

1 Format ( 1 Mars-GRAM Batch version 3.34 - November 1, 1995') SETU 4 

C SETU 5 

***** Mars-GRAM 3 . 4XSETUP . FOR 

C SETU 3 

1 Format (’ Mars-GRAM Batch version 3.4 - April 1, 1996') SETU 4 

C SETU 5 

* * ★ * * 


The climate adjustment factors (CF) and the thermospheric adjustment parameters (DZF = 
deltaZF, dTF = deltaTF, and dTEX = deltaTEX) are added to the common DATACOM. 


***** Mars-GRAM 3 . 34\SETUP . FOR 

& logscale, dustlat , dustlon, dusthgt , radmax, Rref , modpert , alsO , SETU 19 

& intens , iuO , iup, maxf iles SETU 20 

C SETU 2 1 

***** Mars-GRAM 3 . 4 \ SETUP . FOR 

& logscale, dustlat, dustlon, dusthgt, radmax, Rref , modpert , alsO , SETU 19 

Sc intens,iu0,iup,maxfiles,CF(0:5) ,dZF, dTF, dTEX SETU 20 

C SETU 21 

** *** 


Default values are established for all of the input parameters. The default values are those 
that are applicable to the Viking 1 lander site date and time, at the Viking 1 lander surface position. 
With the default values set, only those parameters that differ from their default values need to be 
input in the NAMELIST file. 


***** Mars-GRAM 3 . 34\SETUP . FOR 

DATA IDAY/0 # 31, 59, 90, 120, 151, 181, 212,243, 273,304, 334/ SETU 42 

C SETU 43 

***** Mars-GRAM 3 . 4 \ SETUP . FOR 

DATA IDAY/0 , 31, 59 , 90 , 120 , 151 , 181 , 212 , 243 , 273 , 304 , 334/ SETU 42 

C.. . Establish default values for input parameters SETU 42a 

Data LSTFL, OUTFL/ ' LIST ' , ’OUTPUT 1 / SETU 42b 

C... Default time = Viking 1 Lander SETU 42c 

Data Month, Mday , Myear , Ihr , Imin, Sec /7, 20,76, 12, 30, 0.0/ SETU 42d 

C. . . Default Solar Flux parameters SETU 42e 

Data F107 , Stdl/68 . ,0.0/ SETU 42f 

C... Default position = Viking 1 Lander Site, height = 0 SETU 42g 

Data Flat , Flon, Fhgt , Delhgt , Del la t, Dellon, Deltime/22 . , 48 . ,-0.5, SETU 42h 

& 10. ,3*0./ SETU 42i 

C... CF climate factors default to 1.0 SETU 42 j 

Data CFO , CF5 , CF15 , CF30, CF50 , CF75/6*1 . 0/ SETU 42k 

C... Terms deltaZF, deltaTF and deltaTEX default to 0.0 SETU 421 

Data deltaZF, deltaTF, deltaTEX/3*0 . 0/ SETU 42m 

C... Default random number seed SETU 42n 

Data NR1/1001/ SETU 42o 

C SETU 43 


The climate adjustment factors (CF) and the thermospheric adjustment parameters 
(deltaZF, deltaTF, and deltaTEX) are added to the NAMELIST definition. 
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***** Mars-GRAM 3 . 34\SETUP . FOR 

& NVARX , NVARY , LOG SCALE , FLAT , FLON , FHGT , DELHGT , DELLAT , DELLON , SETU 47 

& DELTIME SETU 48 

SETU 49 

***** Mars-GRAM 3 . 4 \ SETUP . FOR 

& NVARX , NVARY , LOGSCALE , FLAT , FLON , FHGT , DELHGT , DELLAT , DELLON , SETU 47 

& DELTIME / CF0 / CF5,CFl5,CF30,CF50,CF75,deltaZF,deltaTF / deltaTEX SETU 48 

SETU 49 

★ * * ★ ★ 

The climate adjustment factors and thermospheric adjustment parameters are loaded into 
the common DATACOM. 

***** Mars-GRAM 3 . 34\SETUP .FOR 

Read (8, INPUT) SETU 53 

SETU 54 

***** Mars-GRAM 3 . 4XSETUP . FOR 

Read (8, INPUT) SETU 53 

CF ( 0 ) = CFO SETU 53a 

CF ( 1 ) = CF5 SETU 53b 

CF ( 2 ) = CF15 SETU 53c 

CFO) = CF30 SETU 53d 

CF (4) = CF50 SETU 53e 

CF ( 5 ) = CF75 SETU 53f 

Do 5 i = 0,5 SETU 53g 

If (CFO) .le. 0.0) Stop 1 Bad CF value' SETU 53h 

5 Continue SETU 53 i 

dzF = deltaZF SETU 53 j 

dTF = deltaTF SETU 53k 

dTEX = deltaTEX SETU 531 

SETU 54 

***** 
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Appendix B - Reference LIST Output for Mars-GRAM 3.4 


Mars-GRAM Batch version 3.4 - April 1, 1996 
Date = 7/20/1976 Julian Date = 2442980.0 GMT Time =12:30: .0 

F10. 7 flux = 68.0 (1 AU) 25.0 (Mars), standard deviation = .0 

Perturbation model = 3 Starting random number = 1001 

Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = -.50 km ( .00 km) Scale Height = 12.50 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 243.8 K Pressure = 7.479E+02 N/ro**2« 

Density (Low, Avg., High) = 1.500E-02 1.604E-02 1.709E-02 kg/m**3 

Departure, COSPAR NH Mean = -7.6% -1.2% 5.3% 

Density perturbation = 5.72 % of mean value 

Eastward Wind = 4.3 m/s Northward Wind = -.4 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 9.50 km ( 10.00 km) Scale Height = 10.40 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 202.9 K Pressure = 3.058E+02 N/m**2 

Density (Low, Avg., High) = 7.621E-03 7.886E-03 8.154E-03 kg/m**3 

Departure, COSPAR NH Mean = ' 12.9 % 16.8 % 20.8 % 

Density perturbation = -1.39 % of mean value 

Eastward Wind = 1.4 m/s Northward Wind = -5.3 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 19.50 km ( 20.00 km) Scale Height = 9.50 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 185.2 K Pressure = 1.124E+02 N/m**2 

Density (Low, Avg., High) = 3.041E-03 3.175E-03 3.312E-03 kg/m**3 

Departure, COSPAR NH Mean = 10.3 % 15.1 % 20.1 % 

Density perturbation = -1.42 % of mean value 

Eastward Wind = -1.7 m/s Northward Wind = -8.7 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 29.50 km ( 30.00 km) Scale Height = 8.73 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 170.4 K Pressure = 3.812E+01 N/m**2 

Density (Low, Avg., High) = 1.098E-03 1.170E-03 1.244E-03 kg/m**3 

Departure, COSPAR NH Mean = 6.2% 13.2% 20.3% 

Density perturbation = -.33 % of mean value 

Eastward Wind = -4.6 m/s Northward Wind = -13.5 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 39.50 km ( 40.00 km) Scale Height = 8.20 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 160.0 K Pressure = 1.198E+01 N/m**2 

Density (Low, Avg., High) = 3.667E-04 3.915E-04 4.173E-04 kg/m**3 

Departure, COSPAR NH Mean = 2.1% 9.1% 16.2% 

Density perturbation = -.34 % of mean value 

Eastward Wind = -8.0 m/s Northward Wind = -19.7 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 49.50 km ( 50.00 km) Scale Height = 7.68 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 
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Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 149.8 K Pressure = 3.511E+00 N/m**2 

Density {Low, Avg., High) = 1.108E-04 1.226E-04 1.347E-04 kg/m**3 

Departure, COSPAR NH Mean = -3.5% 6.8% 17.3% 

Density perturbation = 5.25 % of mean value 

Eastward Wind = -11.7 m/s Northward Wind = -26.4 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 59.50 km ( 60.00 km) Scale Height = 7.40 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 144.3 K Pressure = 9.703E-01 N/m**2 

Density {Low, Avg., High) = 3.040E-05 3.517E-05 4.002E-05 kg/m**3 

Departure, COSPAR NH Mean = -10.3% 3.8% 18.1% 

Density perturbation = -8.82 % of mean value 

Eastward Wind - -15.9 m/s Northward Wind = -33.7 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height * 69.50 km ( 70.00 km) Scale Height = 7.14 km 

Latitude » 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Temperature = 139.2 K Pressure = 2.579E-01 N/m**2 

Density {Low, Avg., High) = 8.376E-06 9.690E-06 1.103E-05 kg/m**3 

Departure, COSPAR NH Mean = -10.2 % 3.9 % 18.2 % 

Density perturbation = -2.23 % of mean value 

Eastward Wind = -20.4 m/s Northward Wind = -3 9.8 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 79.50 km ( 80.00 km) Scale Height = 7.25 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 128.2 K Zbase * 117.1 km 

Temperature = 134.9 K Pressure = 6.050E-02 N/m**2 

Density (Low, Avg., High) = 1.838E-06 2.343E-06 2.892E-06 kg/m**3 

Departure, COSPAR NH Mean = -24.9 % -4.3 % 18.1 % 

Density perturbation = 3 . 56 % of mean value 

Eastward Wind = -25.4 m/s Northward Wind = -43.8 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 89,50 km ( 90.00 km) Scale Height = 7.16 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase - 128.2 K Zbase = 117.1 km 

Temperature = 132.1 K Pressure = 1.182E-02 N/m**2 

Density (Low, Avg., High) = 3.388E-07 4.663E-07 6.134E— 07 kg/m**3 

Departure, COSPAR NH Mean = -47.3 % -27.4 % -4.5 % 

Density perturbation = -34.81 % of mean value 

Eastward Wind = -38.3 m/s Northward Wind = -45.7 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 99.50 km ( 100.00 km) Scale Height = 7.11 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 128.2 K Zbase = 117.1 km 

Temperature = 129.9 K Pressure = 2.279E-03 N/m**2 

Density {Low, Avg., High) = 6.544E-08 9.113E-08 1.211E-07 kg/m**3 

Departure, COSPAR NH Mean = -61.5 % -46.4 % -28.8 % 

Density perturbation = -9.68 % of mean value 

Eastward Wind = -62.1 m/s Northward Wind = -47.3 m/ s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 109.50 km { 110.00 km) Scale Height = 7 . 11 km 
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Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg. 

Sun Longitude = 108.77 deg. 

Exospheric Temp. = 210.9 K 

Temperature = 128.9 K 

Density (Low, Avg., High) = 
Departure, COSPAR NH Mean = 
Density perturbation = 
Eastward Wind = -103.2 m/s 


West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117.1 km 

Pressure = 4.402E-04 N/m**2 

1 . 269E-08 1 . 768E-08 2.350E-08 kg/m**3 

-71.3 % -60.0 % -46.8 % 

24.11 % of mean value 

Northward Wind = -50.0 m/s 


Time (rel. to TO) = 

Height = 119.50 km { 120.00 
Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg. 

Sun Longitude = 108.77 deg. 

Exospheric Temp. = 210.9 K 

Temperature = 144.9 K 

Density (Low, Avg., High) = 
Departure, COSPAR NH Mean = 
Density perturbation = 
Eastward Wind = -139.0 m/s 


0 sec. ( .000 sols) Ls = 97.0 deg. 

km) Scale Height = 8.09 km 

West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117.1 km 

Pressure = 9.240E-05 N/m**2 

2 . 358E-09 3 . 285E-09 4.366E-09 kg/m**3 

-81.4 % -74.1 % -65.5 % 

-23.70 % of mean value 

Northward Wind = -49.7 m/s 


Time (rel. to TO) = 

Height = 129.50 km ( 130.00 
Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg. 

Sun Longitude = 108.77 deg. 

Exospheric Temp. = 210.9 K 

Temperature = 185.2 K 

Density (Low, Avg., High) = 
Departure, COSPAR NH Mean = 
Density perturbation = 
Eastward Wind = -177.2 m/s 


0 sec. { .000 sols) Ls = 97.0 deg. 

km) Scale Height = 10.57 km 

West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117.1 km 

Pressure = 3.201E-05 N/m**2 

6.286E-10 8 . 759E-10 1.164E-09 kg/m**3 

-84.2 % -78.0 % -70.7 % 

.46 % of mean value 
Northward Wind = -72.5 m/s 


Time (rel. to TO) = 

Height = 139.50 km ( 140.00 
Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg. 

Sun Longitude = 108.77 deg. 

Exospheric Temp. = 210.9 K 

Temperature = 200.9 K 

Density (Low, Avg., High) = 
Departure, COSPAR NH Mean = 
Density perturbation = 
Eastward Wind = -193.9 m/s 


0 sec. { 


.000 sols) 


Ls = 


97.0 deg. 


km) Scale Height = 11.78 km 

West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117,1 km 

Pressure = 1.313E-05 N/m**2 

2 . 3 25E-10 3 . 240E-10 4.306E-10 kg/m**3 

-79.8 % -71.8 % -62.5 % 

7.64 % of mean value 
Northward Wind = -108.6 m/s 


Time (rel. to TO) = 

Height = 149.50 km ( 150.00 
Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg 

Sun Longitude = 108.77 deg 

Exospheric Temp. = 210.9 K 

Temperature = 206.9 K 

Density (Low, Avg., High) = 
Departure, COSPAR NH Mean = 
Density perturbation = 
Eastward Wind = -178.4 m/s 


. 0 sec. ( 


.000 sols) 


Ls = 


97.0 deg. 


km) Scale Height = 12.59 km 

West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117.1 km 

Pressure = 5.780E-06 N/m**2 

9 . 633E-11 1 . 342E-10 1.784E-10 kg/m**3 

-80.4 % -72.7 % -63.7 % 

24.11 % of mean value 

Northward Wind = -157.3 m/s 


Time (rel. to TO) = .0 sec. 

Height = 159.50 km ( 160.00 km) 
Latitude = 22.000 degrees 

Sun Latitude = 25.00 deg. 

Sun Longitude = 108.77 deg. 

Exospheric Temp. = 210.9 K 

Temperature = 209.3 K 

Density (Low, Avg., High) = 4 

Departure, COSPAR NH Mean = 


( .000 sols) Ls = 97.0 deg. 

Scale Height = 13.39 km 

West Longitude = 48.000 degrees 

Mars Orbital Radius = 1.649 AU 
Local Time = 16.05 Mars hours 

Tbase = 128.2 K Zbase = 117.1 km 

Pressure = 2.674E-06 N/m**2 

. 2 17E-11 5 . 875E-11 7.808E-11 kg/m* 

-83.2 % -76.6 % -68.9 % 
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Density perturbation = -46.41 % of mean value 

Eastward Wind = -112.9 m/s Northward Wind = -191.6 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 169.50 km ( 170.00 km) Scale Height = 14.36 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 12 8.2 K Zbase = 117.1 km 

Teirperature = 210.2 K Pressure = 1.298E-06 N/m**2 

Density (Low, Avg., High) = 1.919E-11 2.674E-11 3.554E-11 kg/m**3 

Departure, COSPAR NH Mean = -86.2 % -80.8 % -74.4 % 

Density perturbation = -14.25 % of mean value 

Eastward Wind = -33.4 m/s Northward Wind = -168.4 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 179.50 km ( 180.00 km) Scale Height = 15.64 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 128.2 K Zbase = 117.1 km 

Temperature = 210.6 K Pressure = 6.652E-07 N/m**2 

Density (Low, Avg., High) = 9.081E-12 1.265E-11 1.682E-11 kg/m**3 

Departure, COSPAR NH Mean = -88.8 % -84.4 % -79.3 % 

Density perturbation = -.62 % of mean value 

Eastward Wind = 6.0 m/s Northward Wind = -116.9 m/s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 189.50 km ( 190.00 km) Scale Height = 17.30 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 128.2 K Zbase = 117.1 km 

Temperature = 210.8 K Pressure = 3.617E-07 N/m**2 

Density (Low, Avg., High) = 4.489E-12 6.254E-12 8.312E-12 kg/m**3 

Departure, COSPAR NH Mean = —90.8% -87.1% -82.9% 

Density perturbation = .27 % of mean value 

Eastward Wind = 15.7 m/s Northward Wind = -77.2 m/ s 


Time (rel. to TO) = .0 sec. ( .000 sols) Ls = 97.0 deg. 

Height = 199.50 km ( 200.00 km) Scale Height = 19.35 km 

Latitude = 22.000 degrees West Longitude = 48.000 degrees 

Sun Latitude = 25.00 deg. Mars Orbital Radius = 1.649 AU 

Sun Longitude = 108.77 deg. Local Time = 16.05 Mars hours 

Exospheric Temp. = 210.9 K Tbase = 128.2 K Zbase = 117.1 km 

Temperature = 210.8 K Pressure = 2.092E-07 N/m**2 

Density (Low, Avg., High) = 2.334E-12 3.252E-12 4.322E-12 kg/m**3 

Departure, COSPAR NH Mean = -92.2% -89.1% -85.6% 

Density perturbation = -8.87 % of mean value 

Eastward Wind = 15.5 m/s Northward Wind = -52.9 m/s 
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